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Abstract 

In this paper, we study the inverse scattering of massive charged Dirac fields in the exterior region 
of (de Sitter)-Reissner-Nordstr6m black holes. First we obtain a precise high-energy asymptotic 
expansion of the diagonal elements of the scattering matrix (i.e. of the transmission coefficients) and 
we show that the leading terms of this expansion allows to recover uniquely the mass, the charge and 
the cosmological constant of the black hole. Second, in the case of nonzero cosmological constant, 
we show that the knowledge of the reffection coefficients of the scattering matrix on any interval of 
energy also permits to recover uniquely these parameters. 

1 Introduction 

This paper deals with inverse scattering problems in black hole spacetimes and is a continuation of our 
previous work [31 . Here we shall study the inverse scattering of massive charged Dirac fields that propagate 
in the outer region of (de Sitter)-Reissner-Nordstrom black holes, an important family of spherically 
symmetric, charged exact solutions of the Einstein equations that will be thoroughly described in Section 
[H These spacetimes are completely characterized by three parameters: the mass M > and the electric 
charge Q gM. of the black hole and the cosmological constant A > of the universe. In what follows, these 
parameters will be considered as the "unknowns" of our inverse problem. In fact, the inverse scattering 
problem we have in mind is of the following type: we assume that we are observers living in the exterior 
region of a (dS)-RN black hole, that is the region between the exterior event horizon of the black hole 
and the cosmological horizon when A > 0, or the region lying beyond the exterior event horizon of the 
black hole when A = 0. The geometry of the spacetime in which these observers live is thus fixed in 
some sense. But, what we don't assume however is that these observers know the exact values of the 
parameters M,Q and A "a priori". Hence the natural question we adress is: do such observers have any 
means to measure or characterize uniquely these parameters by an inverse scattering experiment? 

Let us first describe more precisely the exact inverse scattering experiment studied in this paper. A 
direct scattering theory for massive charged Dirac fields has been established in [3] for RN black holes and 
more generally in [l8j for dS-RN black holes. As shown in these papers, the outer region of (dS)-RN black 
holes offers an original geometrical situation from the point of view of scattering theory. These spacetimes 
possess indeed two distinct asymptotic regions, namely either the exterior event horizon of the black hole 
and the cosmological horizon when A > 0, or the event horizon of the black hole and spacelike infinity 
when A = 0, which may have very different geometrical structures. The first consequences concerning the 
propagation properties of Dirac fields are given in the following important result obtained in [3l [18] : the 
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energy of Dirac fields contained in any compact set between the two asymptotic regions tends to zero when 
the time tends to infinity. Therefore, Dirac fields scatter toward these asymptotic regions at late times 
and moreover, they are shown to obey there simple but different equations. From the mathematical point 
of view, two distinct wave operators must be introduced according to the asymptotic region we consider. 
Let us denote for the moment the wave operators corresponding to the part of Dirac fields which scatters 
toward the event horizon of the black hole by W^^.^ and the wave operators corresponding to the part 

of Dirac fields which scatters toward the cosmological horizon or spatial infinity by W^^y These wave 
operators will be precisely defined in Section [2l Now the main result obtained in [3l [18] shows that the 
global wave operators defined by 

w'^ = wf' , + (1.1) 

(—00) (+00)' ^ -' 

exist and are asymptotically complete. This permits to define a global scattering operator S by the usual 
formula 

s = {w+yw-. 

The scattering operator S will be the main object of study of this paper. In fact, we rephrase and 
precise our initial problem in the following way. We assume that our observers have access experimentally 
to the scattering operator S. In particular, we assume that they may measure the expectation values 
of 5, i.e. they can measure any quantities of the form < Sip, (I) > where < .,. > denotes the scalar 
product of the energy Hilbert space Ti. on which S acts and ijj, (j) are any element of Ti.. The question we 
adress is now; is the knowledge of S and any of its related quantities a sufficient information to uniquely 
characterize the parameters M, Q and A of (dS)-RN black holes? 

We can in fact be more precise in the statement of the problem if we remark that the scattering 
operator 5* can be decomposed using (|1.1|1 as 

S = Tl+Tr + L + R, 

where 
and 

R = (W^(+oo))*W^r+oo)> L = 
Each of the terms in S corresponds to a different inverse scattering experiment. For instance, the first 
two terms Tr and (in fact the diagonal elements of S) can be understood as transmission operators. 
Precisely, they correspond to the following experiment: a signal is emitted in the remote past from one 
asymptotic region and is captured in the late future in the other asymptotic region. These two terms 
measure thus the part of a signal which is transmitted from one asymptotic region to the other in a 
scattering process. Conversely, the last two terms L and R (the anti-diagonal elements of S) can be 
understood as reflection operators and correspond to the opposite experiment: a signal is emitted in the 
remote past from one asymptotic region and is captured in the late future in the same asymptotic region. 
Hence these two terms measure the part of the signal which is reflected from an asymptotic region to itself 
in a scattering process. Depending on the point of view of our observers, the quantities of interest will 
be thus the expectation values < TriP, >, < TliP, > and < Lip, (p >,< Rip, (p > oi the transmission 
and reflection operators respectively. In this paper, we shall study in fact two types of inverse problems. 
Firstly, in the two cases of RN black holes (A = 0) and dS-RN black holes (A > 0), we shall prove that 
the parameters M, Q, A are uniquely determined if we assume that the high energies of the transmission 
operators Tr or Tl are known. Secondly, in the case of dS-RN black holes only (A > 0), we shall prove 
the same uniqueness result under the assumption that the reflection operators L or R arc known on any 
(possibly small) interval of energy. 
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Let us now recall the results of [4] where the first kind of inverse problem was adressed in the case of 
Reissner- Nordstrom black holes {i.e. with only the two parameters M, Q unknown and the cosmological 
constant A equal to 0). Using the direct scattering theory for massless Dirac fields obtained in [51 
and a high energy asymptotic expansion of the expectation values < Tji^, </) > or < T^ip, 4> > (^-s defined 
above), a partial answer was then given: the mass M and the modulus of the charge \Q\ are uniquely 
determined from the leading terms of this high energy asymptotic expansion. Note that the indecision of 
the sign of the charge is not surprising in that case since the propagation of massless Dirac fields is only 
influenced by the geometry of the black hole which in turn only depends on \Q\ (see the expression of 
the metric (I2.2p in Section [J). Moreover, it was mentioned in [1] (see also [H] where a similar problem 
was studied) that the same method couldn't be applied to uniquely recover the parameters from the 
high energies of the reflection operators R oi L. The relevant quantities are in that case indeed non 
measurable. In this paper we continue our investigation and improve our results in several directions. 

In Section [3l we reconsider the case A = corresponding to RN black holes but study the inverse 
scattering of massive charged Dirac fields instead of massless Dirac fields. Using the same approach than 
in [4], we show that the mass M as well as the charge Q are uniquely determined by the leading terms 
of the high energy asymptotic expansion of the transmission operators Tr or T^. In fact, the advantage 
of considering massive charged Dirac fields is that an explicit term associated to the interaction between 
the electric charge of the fields and that of the black hole appears in the equation and allows to recover 
Q and not \Q\. From the mathematical side, the analysis turns out to be much more involved than in 
[4] because of two main reasons. First, massive Dirac fields have completely distinct behaviours when 
approaching the different asymptotic regions. At the event horizon of the black hole for instance, the 
attraction exerced by the black hole is so strong that massive Dirac fields behave as massless Dirac fields. 
The dynamic there is very simple and will be shown to obey a system of transport equations along the 
null radial geodesies of the black hole. This is a consequence of the particular geometry (of hyperbolic 
type) near the event horizon (and more generally near any horizons). Conversely, RN black holes are 
asymptotically flat at spacelike infinity. There, the fields simply behave like massive Dirac fields in the 
usual Minkowski spacetime and the mass of the fields, slowing down the propagation, plays an important 
role. In consequence the dynamics near the two asymptotic regions are quite different and must be 
treated separatly. The second kind of difficulty comes from the appearance of long-range terms in the 
equation but only in one asymptotic region: spacelike infinity. This entails new technical difficulties such 
as a modification of the standard wave operators at infinity and we need to work harder to obtain the 
high energy asymptotic expansion of the transmission operators. We want to emphasize at last that the 
model studied in this part can be viewed as a good intermediate model before studying the same inverse 
problem in the more complicated geometrical setting of Kerr black holes. As shown in [13] indeed, the 
appearance of long-range terms in the equation (even for massless Dirac fields) is compulsory in that case 
as a side effect of the rotation of the spacetime. 

In Section [4l we consider the case of nonzero cosmological constant A > 0, that is de Sitter-Reissner- 
Nordstrom black holes and we have three unknown parameters M,Q,A a priori. The two asymptotic 
regions are the event horizon of the black hole and the cosmological horizon. Near these regions, massive 
Dirac fields behave as massless Dirac fields and as before, their propagation obeys essentially a system 
of transport equations along the null radial geodesies of the black hole. However, different oscillations 
appear in the dynamics near these two horizons, once again due to the interaction between the charge 
of the field and that of the black hole. In consequence, Dirac fields evolve according to slightly different 
dynamics in that case too. In Subsection [131 using the results of the previous part, we shall obtain a high 
energy asymptotic expansion of the transmission operators Tr and Tl and again, we shall prove that the 
parameters M, Q and A are uniquely characterized by the leading terms of this asymptotic expansion. 
Then we consider inverse scattering experiments based on the knowledge of the reflection operators R 
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or L on a (small) interval of energy. As already mentioned, a high energy aymptotic expansion of these 
reflection operators doesn't give any information and can't be used to solve the inverse problem. To 
study this case, we follow instead the usual stationary approach of inverse scattering theory on the line. 
We refer for instance to the review by Faddeev [8] and to the important paper by Deift and Trubowitz 
[6] for a presentation of the method for Schrodinger operators and to the nice paper [T] for a recent 
application to Dirac operators (see also [El HI])- In Subsection 14.21 we shall first obtain a stationary 
representation of the scattering operator S in terms of the usual transmission and reflection " coefficients" 
(note that these turn out to be matrices in our case). This is done after a serie of simplifications of 
our model which happens finally to reduce to a particular case of the model studied in [T]. Then we 
use the analysis of [1] , namely a classical Marchenko method based on a carefuU study of the stationary 
solutions of the corresponding Dirac equation, to prove the following result: the knowledge of one of the 
reflection operators L or i? at all energies is enough to uniquely characterize the parameters M, Q and 
A. Eventually, we improve this result seeing that, in our model, the reflection operators R or L are in 
fact analytic in the energy variable on a small strip containing the real axis. Hence it is enough to know 
i? or L on any interval of energy in order to uniquely know them for all energies. Applying the result of 
[1], this leads to the uniqueness of the parameters in that case too. 

We finish this introduction saying a few words on the main technical tool used in Sections [3] and |4] 
to prove our uniqueness results from the high energies of the transmission operators Tr ot Tl- These 
are based on a high-energy expansion of the scattering operator S following an approach introduced by 
Enss and Weder in [7] in the case of multidimensional Schrodinger operators. (Note that the case of 
multidimensional Dirac operators in flat spacetime was treated later by Jung in |17|). Their result can 
be summarized as follows. Using purely time-dependent methods, they showed roughly speaking that 
the first term of the high-energy expansion of S is exactly the Radon transform of the potential they 
are looking for. Since they work in dimension greater than two, this Radon transform can be inversed 
and the potential thus uniquely recovered. In our problem however, due to the spherical symmetry of 
the black hole, we are led to study a family of one dimensional Dirac equations and the above Radon 
transform simply becomes an integral of a one-dimensional function, hence a number, and cannot be 
inversed. Fortunately in our models, it turns out that this integral can be explicitely computed and 
gives in general already a physically relevant information. Nevertheless, it is not enough to uniquely 
characterize all the parameters of the black hole. In fact, we need to calculate several terms of the 
asymptotic (and thus obtain several integrals) to prove our result. To do this, we follow the stationary 
technique introduced by one of us [21] which is close in spirit to the Isozaki-Kitada method used in 
long-range scattering theory [16] . The basic idea is to replace the wave operators (and thus the scattering 
operator) by explicit Fourier Integral Operators, called modifiers, from which we are able to compute 
the high-energy expansion readily. The construction of these modifiers and the precise determination of 
their phases and amplitudes will be given in a self-contained manner in Section [3J Note also that the 
similar results proved in our previous paper [4] couldn't be applied directly to our new model because 
of the presence of long-range terms in the equation. At last we mention that, while this method was 
well-known for Schrodinger operators and apphed successfully to various situations (see [21211 [22l|23] ) , it 
has required some substantial modifications when applied to Dirac operators, essentially because of the 
matrix-valued nature of the equation. To deal with these difficulties, we made an extensive use of the 
paper by Gatel and Yafaev [9] where a direct scattering theory of massive Dirac fields in flat spacetime 
was studied and modiflers were constructed. 
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2 (De Sitter)-Reissner Nordstrom black holes and Dirac equa- 
tion 



In this section, we describe the geometry of the exterior regions of (de Sitter)-Reissner-Nordstr6m black 
holes. In particular, we emphasize the different properties of the asymptotic regions mentioned in the 
introduction, clearly distinguishing between the cases of zero and nonzero cosmological constant A. We 
then express in a synthetic manner the equations that govern the evolution of massive charged Dirac 
fields in these spacetimes. We end up this section recalling the known direct scattering results of [3l [18] 
and introducing the scattering operator S. 

2.1 (De Sitter)-Reissner-Nordstr6m black holes 

In Schwarzschild coordinates a (de Sitter)-Reissner-Nordstr6m black hole is described by a four dimen- 
sional smooth manifold 

X =Mt x]R+ X Si, 

equipped with the lorentzian metric 

g ^ F(r)df - F{r)-^dr^ -r^duj^, (2.1) 

where 



= — + ^-— , (2.2) 



2Af Q2 ^^2 
r 3 

and duj^ — dO^ + sir? 9 dip^ is the euclidean metric on the sphere S^. The constants M > 0, Q G M 
appearing in (|2.2p are interpreted as the mass and the electric charge of the black hole and A > is 
the cosmological constant of the universe. Observe that the function (|2.2|) and thus the metric (|2.1[) do 
not depend on the angular variables 9,(p E S"^ reflecting the fact that dS-RN black holes are spherically 
symmetric spacetimes. 

The family {A4 , g) are in fact exact solutions of the Einstein-Maxwell equations 



Gf_iu — SnTf^u, Gf_,i, — Rfii, + -Rg^i, + Ag^jy, (2.3) 



Here Gfii,,Rp,i, and R denote respectively the Einstein tensor, the Ricci tensor and the scalar curvature 
of {A4 , g) while Tp,y is the energy- momentum tensor 

T^u = ^ {F^pF/ - -^g^^Fp^Ff"^) , (2.4) 

where F^^, is the electromagnetic two-form solution of the Maxwell equations V^F^p — 0, V[^i^yp] = 
and given here in terms of a global electromagnetic vector potential 

Fpu = V[pA,], A.dx" = -'^dt. (2.5) 



We point out that any spherically symmetric solutions of the Einstein equations (|2.3p - ()2.5|l must belong 
(at least locally) to the family of dS-RN black holes defined by (12. ip and ()2.2p . This is a well-known 
uniqueness result due to Birkhoff (see for instance [14j). In particular, the results contained in this paper 
apply to this extended class of spacetimes. 
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The metric g has two types of singularities. Firstly, the point {r = 0} for which the function F is 
singular. This is a true singularity or curvature singularity q Secondly, the spheres whose radii are 
the roots of F (note that the coefficient of the metric g involving F^^ explodes in this case). We must 
distinguish here two cases. When the cosmological constant is positive A > and small enough, there 
are three positive roots < ^- < < < +oo . The spheres of radius r_ , rg and r^ are called 
respectively Cauchy, event and cosmological horizons of the dS-RN black hole. When A = 0, the number 
of these roots depends on the respective values of the constants M and Q. In this paper we only consider 
the case M > \Q\ for which the function F has two zeros at the values r_ = M — — and 

ro — M + ^ AP — The spheres of radius and tq are called respectively the Cauchy and event 
horizons of the RN black hole. In both situations, the horizons are not true singularities in the sense 
given for {r = 0}, but in fact coordinate singularities. It turns out that, using appropriate coordinate 
systems, these horizons can be understood as regular null hypersurfaces that can be crossed one way but 
would require speeds greater than that of light to be crossed the other way. We refer to [13] and for 
a introduction to black hole spacetimes and their general properties. 

As mentioned in the introduction, we shall consider in this paper inverse scattering experiments made 
by observers living in the exterior region of a (dS)-RN black hole, that is the region {rp < r < r+} 
when A > or the region {rp < r < +00} when A = 0. It is thus important to understand the roles 
of the horizons as the natural boundaries of the exterior region. In Schwarzschild coordinates, it turns 
out that they are asymptotic regions of spacetime. Precisely, this means that they are never reached in 
a finite time t by incoming and outgoing null radial geodesies, i.e the trajectories followed by classical 
light-rays aimed radially at the black hole and either at the cosmological horizon if A > or at infinity if 
A = 0. To see this point more easily, we introduce a new radial coordinate x, called the Regge- Wheeler 
coordinate, which has the property of straightening the null radial geodesies and will, at the same time, 
greatly simplify the later analysis. Observing that for all A > the function F{r) in the metric 
remains always positive in the exterior region, it can be defined implicitely by the relation 



dr 
dx 



^ Fir) > 0, (2.6) 



or explicitely, by 

1 

where the quantity 



log(r - ro) - r - p-)dy] + C, (2.7) 
= \F\ro) > 0, 

is called the surface gravity of the event horizon and C is any constant of integration. Note that, when 
A > 0, the Regge- Wheeler variable could be also defined explicitely by 

log(r+ - r) - r + ^My] + C, (2.8) 
Jr r+-y F{y)' J 



1 

2k] 



where the quantity 

is called the surfave gravity of the cosmological horizon. Moreover, in the case A = 0, the expression 
(|2.7p simplifies as 

1 

x^r+— \og{r - ro) H log(r - r_) + C. (2.9) 

2ko ro - r_ 



^It means that certain scalars obtained by contracting the Riemann tensor blow up when r — > 0. 
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In the coordinate system {t,x,uj), it is easy to see from the logarithm in ()2.7p and (|2.9p and the 
positive sign of kq that the event horizon {r — vq} is pushed away to {x — —00} for aU A > 0. Similarly 
it follows from (|2.8p and the negative sign of k+ that the cosmological horizon {r — r+} is pushed away 
to {a; = +00} when A > 0. Hence in any case the Regge- Wheeler variable x runs over the full real line 
M. Moreover, by (|2.6|) . the metric takes now the nice form 

g = F{r){dt^ - dx^)-r^du;'^, (2.10) 

from which it is immediate to see that the incoming and outgoing null radial geodesies are generated by 
the vector fields ^ ± ^ and take the simple form 

j^it)^{t,XQ±t,LOo), teR, (2.11) 

where {xq,ujo) E M. x S'^ are fixed. These are simply straight lines with velocity ±1 mimicking, at least 
in the t — x plane, the situation of a one-dimensional Minkowski spacetime. At last, using (|2.1ip . we can 
check directly that the event horizon and the cosmological horizon (when A > 0) are asymptotic regions 
of spacetime in the sense given above. 

From now on we shall only consider the exterior region of dS-RN black holes and we shall work on 
the manifold S = M( x S with E = R^; x S'^, equipped with the metric (|2.10p . Such a manifold B is 
globally hyperbolic meaning that the foliation Et = {t} x S by the level hypersurfaces of the function t, 
is a foliation of B by Cauchy hypersurfaces (see [26j for a definition of global hyperbolicity and Cauchy 
hypersurfaces). In consequence, we can view the propagation of massive charged Dirac fields as an 
evolution equation in t on the spacelike hypersurface S, that is a cylindrical manifold having two distinct 
ends: {x — —00} corresponding to the event horizon of the black hole and {x — -|-cxd} corresponding to 
the cosmological horizon when A > and to spacelike infinity when A = 0. Note that the geometries 
of these ends are distinct in general. The event and cosmological horizons are indeed exponentially large 
ends of S whereas spacelike infinity is an asymptotically flat end of E (in the latter, observe that the 
metric (|2.2|) tends to the Minkowski metric expressed in spherical coordinates when r +00). The 
difference between these geometries will be easily seen from the distinct asymptotic behaviours of Dirac 
fields near these regions given in the next subsection. 

2.2 Dirac equation and direct scattering results 

Scattering theory for massive charged Dirac fields on the spacetime B has been the object of the papers 
[3l[T8]. We briefiy recall here the main results of these papers. In particular, we use the form of the Dirac 
equation obtained therein. 

First, the evolution equation satisfied by massive charged Dirac fields in B can be written under the 
Hamiltonian form 

idt^ = Hi};, (2.12) 
where ip \s a. 4-components spinor belonging to the Hilbert space 

?^ = L2(KxS'2;C*), 

and the Hamiltonian H is given by 

H = T'^D^ + a{x)Ds2 + b{x)T^ + c{x). (2.13) 

Here we use the following notations. The symbol stands for —idx whereas Dgi denotes the Dirac 
operator on S"^ which, in spherical coordinates, takes the form 

Ds^ = ^^V\^e + ^) - -^r^S^. (2.14) 
2 sm 
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The potentials a, 6, c are scalar smooth functions given in terms of the metric p.ip by 



a{x) 



b{x) = m^/F{r), 



c{x) = , 

r 



(2.15) 



where m and q denote the mass and the electric charge of the fields respectively. Finally, the matrices 
r^, r^, F'^, F° appearing in (|2.13p and (|2.14p are usual 4x4 Dirac matrices that satisfy the anticommutation 
relations 

FT^' + F^F* = 2(5yld, Vi, j = 0, ..,3. (2.16) 

Second, we use the spherical symmetry of the equation to simplify further the expression of the Hamil- 
tonian H. Since, the Dirac operator D52 has compact resolvent, it can be diagonalized into an infinite 
sum of matrix- valued multiplication operators. The eigenfunctions associated to D52 are a generalization 
of the usual spherical harmonics called spin-weighted sphericl harmonics. We refer to I.M. Gel'Fand and 
Z.Y. Sapiro [TU] for a detailed presentation of these generafized spherical harmonics and to [21 US] for an 
application to our model. There exists thus a family of eigenfunctions -F,' of Dg2 with the indexes {l,n) 
running in the setX= {{l,n),l-\^\ eN,l-\n\e N} which fo rms a Hilbert basis of L?{S^; C*) with the 
following property. The Hilbert space 7i can then be decomposed into the infinite direct sum 



(i,n)el 



where Hin = F^(I 



) F,j is identified with (I 



and more important, we obtain the orthogonal 



decomposition for the Hamiltonian H 



H 



H- 

(i,n)el 



In 



with 



= T^D^ + ai{x)T'^ + 6(x)F° + c{x), 



(2.17) 

fin 



and ai{x) ~ ~a{x){l + i). Note that the Dirac operator 1)52 has been replaced in the expression of 
by —{1+ 5)r^ thanks to the good properties of the spin- weighted spherical harmonics F^. The operator 
ij'" is a selfadjoint operator on 7i/„ with domain Z3(i/'") = _ff^(]R;C'*). Finally we use the following 
representation for the Dirac matrices F^, F^ and F° appearing in (|2.17p 



Fi 
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(2.18) 



In this paper it will be often enough to restrict our analysis to a fixed harmonic. To simplify notations 
we shall thus simply write 7i, H and a{x) instead of TLim H'^" and ai{x) respectively and we shall indicate 
in the course of the text whether we work on the global problem or on a fixed harmonic. 

Let us summarize now the direct scattering results obtained in [3l I18j . It is well known that the main 
information of interest in scattering theory concerns the nature of the spectrum of the Hamiltonian H. 
Our first result goes in this sense. Using essentially a Mourre theory (see [H]), it was shown in j3| [T8| 
that, for all A > 0, 

cTppiH) - 0, asi„g(i/) - 0. 

In other words, the spectrum of H is purely absolutely continuous. In consequence, massive charged 
Dirac fields scatter toward the two asymptotic regions at late times and they are expected to obey 
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simpler equations there. This is one of the main information encoded in the notion of wave operators 
that we introduce now. 

We first treat the case A = corresponding to RN black holes. From (I2.2p and (|2.9p . the potentials 
a, 6, c have very different asymptotics as x — > ±cxd (according to our discussion above this reflects the 
fact that the geometries near the two asymptotic regions are very different). At the event horizon, there 
exists a > such that 

|a(x)|, |c(a;)-co| = 0(e"^), x ^ (2.f9) 

where the constant cq is given by (see (|2.f 5p ) 

qQ 

Co = — . 

fo 

Hence, the potentials a,b,c — cq are short-range when x — > —oo and we can choose the asymptotic 
dynamic generated by the Hamiltonian Hq = T^D^ + Cq as the comparison dynamic in this region. The 
Hamiltonian Hq is a selfadjoint operator on Ti. with its spectrum covering the full real line, i.e. (j{Hq) = R. 
Note finally that due to the simple diagonal form of the matrix F^, the comparison dynamic e~'^*^° is 
essentially a system of transport equations along the curves x±t, that is the null radial geodesies of the 
black hole. 

Conversely at infinity, the potentials a, 6, c have the asymptotics 

\a{x)\, \b(x)-ml \c{x)\^0{-), x +oo. (2.20) 

X 

The potentials a,b — m, c are thus long-range potentials having Coulomb decay when x ~f +oo. The 
asymptotic dynamic there is generated by the Hamiltonian i/™ = T^Dx+mT^, a classical one-dimensional 
Dirac Hamiltonian in Minkowski spacetime. The Hamiltonian iJ™ is a selfadjoint operator on Ti. and 
its spectrum has a gap, i.e. a{H^) = (— oo,— m) U {+m, +oo). Contrary to the preceding case, the 
asymptotic dynamic e~^*'^° cannot be used alone as a comparison dynamic because of the long-range 
terms, but must be (Dollard)-modified. 

In order to define this modification and for other use, we need to introduce the classical velocity 
operators 

associated to the Hamiltonians Hq and iJ™ respectively. The classical velocity operators are selfadjoint 
operators on Ti and their spectra are simply a{T^) = {—1, -f 1} and cr(Vm) = [—1, -1-1]. Let us also denote 
by P± and P^^ the projections onto the positive and negative spectrum of F^ and Vm, i.e. 

P± = 1r±(F1), P4- = 1r±(V„). 

As shown in [3], a great interest of these projections is that they permit to separate easily the part of the 
fields that propagate toward the event horizon and the part of the fields that propagate toward infinity. 
They will be used in the definition of the wave operators below. Moreover, the classical velocity operator 
Vm enters in the expression of the DoUard modified comparison dynamic at infinity proposed in [3] and 
given by 

U{t) = e"**-^o"e"*-'o [(f'('*V™)-m)m(ffn")"'+c('*V„)]ds_ (2.21) 

Let us make here two comments. First, the potential a(a;)F^ turns out to be a "false" long-range term. 
This is clear from (I2.2ip where the asymptotic dynamic e~'*^t) has been modified by an extra phase 
which only involves the long-range potentials b and c. We refer to [3] for an explanation of this particular 
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point. Second, we shall propose in Section [3] a new time-independent modification of the comparison 
dynamic e~**^o which will be a direct byproduct of our construction of modifiers in the spirit of Isozaki- 
Kitada's work [TBj. This new modification will be shown to be equivalent to the DoUard modification 
((2:211) in Theorem [331 

We are now in position to introduce the wave operators associated to H. At the event horizon, we 
define 

W^_^^^s- ^lim e'*"e-'*"°P^, (2.22) 



whereas at infinity, we define 

W^. , = ,s - 

t — *±oo 

Finally, the global wave operators are given by 



W^ctoo) = « - . Ii2l„ e"«[/(t)P£'. (2.23) 



W^-Ko.)+Ko.) (2-24) 

Note here our use of the projections P± and P™ to separate the part of the field propagating toward the 
event horizon to the part of the field propagating toward infinity. In fact without these projections, the 
wave operators (|2.22p and (|2.23p wouln't exist at all. More precisely the main result of [3] is 

Theorem 2.1 The wave operators W^^y W^^^ and exist on Ti.. Moreover, the global wave 

operators are partial isometries with initial spaces Ti-^^at — ^t(^) + P±{'H) and final space Ti.. In 
particular, are asymptotically complete, i.e. Ran = Ti. 

As a direct consequence of Theorem 12.11 we can define the scattering operator S by the usual formula 

S^{W+yW~. (2.25) 

It is clear that 5" is a well-defined operator on H and a partial isometry from Tig^at into Tig^at ■ 

We now treat the case A > corresponding to dS-RN black holes wich turns out to be a little bit 
more symmetric at the two (event and cosmological) horizons. According to (|2.2p . (|2.7p and (|2.8p . the 
potentials a, b, c have the following asymptotics as a; — *■ ±oo. There exists a > such that 

|a(x)|, |6(x)| =0(e-"l-l), \x\ ^ ^, (2.26) 

and 

\c{x)-co\ = 0(e°^), x^-oo, (2.27) 

\c{x)-c+\ = 0(e""^), x-^+oo, (2.28) 

where the constants cq and c+ are given by (see (|2.15p ) 

Co = — , c+ = — . (2.29) 

ro r+ 

Hence, the potentials a, b are short-range when x — > ±oo and c — cq and c — c+ are short-range when 
X —oo and x +oo respectively. At the event horizon, we choose as before the asymptotic dynamic 
generated by the Hamiltonian Hq = T^D^ + cq as the comparison dynamic while, at the cosmological 
horizon, we choose the asymptotic dynamic generated by the Hamiltonian iJ+ = T^D^ + c+ as the 
comparison dynamic. The Hamiltonians Hq and i7+ are clearly selfadjoint operators on Ti. and their 
spectra are exactly the real line, i.e. (7{Hq) = a{H^) — R. We observe eventually that the dynamics 
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^-itHa g^j^jj g-it-ff+ g^j.g essentially a system of transport equations along the null radial geodesies of the 
black hole but they differ by the distinct oscillations e^'*'^" and e~'*'^+. 

We need the classical velocity operators associated to and in order to separate the part of 
the fields that propagate toward the event horizon and the part of the fields that propagate toward the 
cosmological horizon, ft turns out that they are equal to Vo — in both cases and the associated 
projections onto the positive and negative spectrum are still P±. Thus we can introduce the wave 
operators as before. At the event horizon, we define 

Wfi ,=s- lim e'*"e-'*"°P^, (2.30) 
and at the cosmological horizon, we define 

K^)-^-,}^^^^''"^-''"^P±- (2-31) 
Finally, the global wave operators are given by 

VF± = VF±_)+M^±^) (2.32) 

The main result of TS] is 

Theorem 2.2 The wave operators W^^y ^^oo) '^^'^ exist on Ti.. Moreover, the global wave 

operators are isometries on H.. In particular, are asymptotically complete, i.e. Ran = Ti. 

Thanks to Theorem 12. 21 we can define the scattering operator S as in (|2.25p by = {W'^)*W^ which is 
a well-defined isometry on Ti. 

We deduce from the previous discussion that, for all A > 0, the scattering operator is a well-defined 
operator on H. For all G Ti, we shall consider in the following the expectation values of S, given 
by < Sij], (j) >, as the known data of our inverse problem. Moreover, using (I2.24p and (|2.32|) . we observe 
that these expectation values can be decomposed into 4 natural components 

<S^,(t)> ^ < W-^, W+<j) > = < Tnij, <j)> + < Tl^, (/)> + < LV', (/)> + < i?VA </> >, 

where 

< Tn^,^>^< Vt^(-^)^, W^(+_)<^ >, < T^^,0 >=< >, (2.33) 

<L7^,0>=<iy(l^)^,W^+^)</'>, <i?V^,0>=<M/(-^)V,W^(Xoo)'^> • (2.34) 

It follows from our definitions of the wave operators (|2.22p . (|2.30p and (|2.23p . (I2.3ip that the previous 
qantities can be interpreted in terms of transmission and reflection between the different asymptotic 
regions, i.e. {x — — oo} for the event horizon of the black hole and {x = -foo} for either spacelike infinity 
if A = 0, or the cosmological horizon if A > 0. For instance, < Taip, (f) > corresponds to the part of a 
signal transmitted from {a; = +00} to {x — —00} in a scattering process whereas the term < TliP, 4> > 
corresponds to the part of a signal transmitted from {a; — —00} to {x ~ +00}. Hence Tr stands for 
"transmitted from the right" and for "transmitted from the left". Conversely, < Lip, (j) > corresponds 
to the part of a signal reflected from {a; = — cx)} to {x = —00} in a scattering process whereas the term 
< Rip, (p > corresponds to the part of a signal reflected from {x = +00} to {x — -l-oo}. 
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3 The inverse problem when A = 



In this section, we study the inverse problem at high energy in the case A = that corresponds to RN 
black holes. Let us recall here that all the results and formulae given hereafter are always obtained 
on a fixed spin-weighted spherical harmonic. Therefore the notations H, H, a{x) are a shorthand for 
TCin, -ff'", ai{x) defined in the preceding Section. In order to state our main result, we make two assump- 
tions. 



Assumption 1 : We assume that our observers may measure the high energies of the transmitted operators 
Tji or Tl. Precisely, we assume that one of the following functions of A G M 

Fi{\) =< TRe'^^^,e'^^<f> >, G/(A) =< Tic'^^V, e*^"0 >, 

are known for all large values of A, for alH € N where I indexes the spin-weighted spherical harmonics 
and for aU TP,(l)Gn with t/., </> e C^(R; C^). 

Assumption 2: We also assume that the mass m and the charge q of the Dirac fields considered in these 
inverse scattering experiments are known and fixed. Moreover we assume that q ^ since the case q — 
is similar to the one treated [4]. 

The main result of this section is now summarized in the following Theorem 

Theorem 3.1 Under assumptions 1 and 2, the parameters M and Q of the RN black hole are uniquely 
determined. 

Following our previous paper [4], the proof of Theorem 13 . 1 1 will be based on a high-energy asymptotic 
expansion of the functions Fi{X) and G;(A) when A +oo. Precisely we shall prove the following 
formulae: 

Theorem 3.2 (Reconstruction formulae) Let ip,(l> E C(j"(M;C^). Then for A large, we obtain 

FiiX) = <e(a;)P_i/^,P_0>+^ <^(a;)P_^,F_0> + 0(A-2), (3.1) 
Gi{X) = <eix)P+^,P+(l)>-^<Aix)F+ij,F+(t>> + 0{\-^), (3.2) 

where 9{x) and A{x) are multiplication operators given by 

A{x)^Q{x)n af{s)ds+ b^{s)ds+ {b{s)-m)^ds + m^xy (3.4) 



Remark 3.1 In Theorem \3.2l we have once again emphasized the dependence of the functions Fi{X) and 
G/(A) on the parameter I since the reconstruction formulae 113. 1\} and iS. 2^) can be derived if we work on 
a fixed spin-weighted spherical harmonic only. Nevertheless, as indicated in Assumption 1 we shall need 
to know these formulae on all spin-weighted spherical harmonics, hence for all I € N, in order to prove 
the uniqueness result stated in Theorem \3.1[ 
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Remark 3.2 In the reconstruction formulae of Thm \3.^ the constant terms j'^^[c{s) — co]ds in iS. 3\) 

and b'^{s)ds + /o^°°(f'(s) — m)'^ds in {S.4\l , that may appear unnatural at first sight since they depend 
explicitely on the particular value of the Regge- Wheeler variable x, are in fact due to our particular 
choice of Dollard modification in the definition of the modified wave operators W^^^ . Recall here indeed 
that there is no canonical choice for the (necessary) modifications entailed by the presence of long-range 
potentials at infinity. This point can be clearly seen for instance from the Isozaki-Kitada modifications 
-constructed in the next Subsection- whose phases are defined only up to a constant of integration (see 
I13.26\) and Remark \3.3\ after it). We emphasize moreover that this constant of integration are not phys- 
ically measurable and we can check indeed that they do not play any role in our proof of the uniqueness 
of the parameters. 

We now explain our strategy to prove Theorem l3.2l Using (|2.22p . (|2.23[) . (|2.33[) and the fact that e^^^ 
corresponds to a translation by A in momentum space, we first rewrite Fi(X) and G/(A) as follows 

Fi{X) = <M/-^)(A)^,W-+^)(A)0>, (3.5) 
Gi{X) = <M/-^)(A)7/;,W^+^)(A)0>, (3.6) 



with 



(+oo)V I ( + oo) t->±oo * 



where we use the notations 

H[X) = V^{D., + A) + a{x)V^ + 6(a;)r° + c{x), H^[X) = V\D., + A) + co, 
Ho"(A) = ri(i?, + A)+mr°, V„(A) = + A) (i7o"(A))"\ P^^^ = 1m± (V„,(A)), 
X(t,A)= / [(6(sV™(A))-m)m(i7^"(A))-i+c(.sV,„(A))]ds. 



In order to obtain an asymptotic expansion of the functions -Fi(A) and G';(A), it is thus enough to obtain 
an asymptotic expansion of the A-shifted wave operators W^^^-^{X). To do this, we follow the procedure 
exposed in [5T1 1^ , procedure inspired by the well-known Isozaki-Kitada method developed in the 
setting of long-range stationary scattering theory. It consists simply in replacing the wave operators 
^(±oo)(^) "well-chosen" energy modifiers J^^^-j(A), defined as Fourier Integral Operators (FIO) with 
explicit phases and amplitudes. Well-chosen here means practically that we look for J^j^^-^ (A) satisfying 
for A large enough 

I^±_)(A)V = ^hm^e"^(^)j±^)(A)e-"^°(^)p^V, (3.7) 
^(too)(A)V' = ^hm^e"^(^)jf+^)(A)e-*«""(^)p™'^V, (3.8) 

and 

ll(W^(loo)W-4oo)(A))^ll =0(A-2), (3.9) 



13 



for any fixed V e H such that -0 e C^{R] C^). Note that the decay 0{X^'^) in jSJ]) could be improved 
to any inverse power decay but turns out to be enough to our purpose here. In particular if we manage 
to construct such J^^.-^ (A) satisfying (|3.9p then we obtain by p.Sp and p.6p 

^;(A) - < J(-^)(A)V,J+^)(A)^> + 0(A-2), (3.10) 

from which we can calculate the first terms of the asymptotics easily. 

Let us here give a simple but useful result which allows us to simplify slightly the expressions of (|3.7p 
and ((3:81) . 

Lemma 3.1 For all f e R*, set 

ly^iO = ±sgn{^)^e + m^- (3.11) 

Then, for all with supp C R* , 

e-'*"^P^ip = e-'^^^^^^^P^i). (3.12) 

Moreover, 

Proof: The Fourier representation of the operator is F^^ + mF° and has precisely one positive 
eigenvalue y^^ + ro^ and one negative eigenvalue — y^^ + rr?. Similarly, the Fourier representation of 
the classical velocity operator V,„ is ^i^^^i i^^S, + rnT'^). Hence, for ^ > 0, P™ is the projection onto 
the positive spectrum of F^^ + mF° and is the projection onto the negative spectrum of F^^ + mF*^. 
For ^ < 0, it is the opposite. This implies immediately p.l2p . Finally the equality (|3.13p is a direct 
consequence of the definitions of Hq and P± . ()■ 

According to Lemma [3.1[ the projections P± and P™ allow us to "scalarize" the Hamiltonians Hq 
and iJ™ in the expressions (|3.7p and p.Sp of W^^JX). Precisely these expressions read now 



^(^oo)(A)V' = ^Hm^e**^(^)j±^)(A)eT^*(^-+^)-**^°P^V, (3.14) 
K+oo)W^ = ^Um^e**^(^)j±^)(A)e-^*''*(^-+^)p£''V- (3.15) 

This minor simplification will be important in the forthcoming construction of the modifiers "^(±00) ("A) 
Before entering into the details, let us give a hint on how to construct the modifiers J(±oo)(^'^ ^ priori 



defined as FIOs with "scalar" phases (yS^^^ (x, ^, A) and "matrix-valued" amplitudes p^^j(x,^, A), i.e. 
defined for all ?A e H by 

J(1^)(A)V ^ -^^e^^f±=^)(-«-%^)(x,C, A)^(Ode 
If we assume for instance that (I3.15P is true then we easily get 

j-±00 

(<+oo)(A) - Jt+oo)im = * e^*^(^)C± ^)(A)e-*'' (^^+^)p£''Vdt, (3.16) 

where 

^^f+oo)(A) := i?(A) J± (A) - J± (A)i.±(7^, + A), (3.17) 
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are also FIOs with phases </?^^^ (x, ^, A) and amplitudes c^^^(x,^,A). From p.l6p we get the simple 
estimate ^ 

mt+oo)W - Jt+oo)WM <jj l|Cf+^)(A) e-*''"(^^+^)pi"'V||di. (3.18) 

In order that (13. 9|) be true it is then clear from (j3.18|) that the FIOs C^^^j(A) have to be "small" in 
some sense. Precisely we shall need that the amplitudes c^^^(a;,^, A) be short-range in the variable 
X at infinity (i.e. when x — > +00) and of order 0(A^^) when A — > +00. Note here the role played 
by the projections P™'^ which allow us to consider the part of the Dirac fields that propagate toward 
infinity. This explains why the amplitudes €^^^(2:,^, A) must short-range in the variable x only at 
infinity. Similarly, for the construction of the modifiers J^_^){^)^ we shall require that the amplitudes 
''f^oo) (^' ^' corresponding operators C^_^-^ (A) be short-range in the variable x only at the event 

horizon {i.e. when x — )■ —00) and of order 0(A^^) when A — > +00. 

3.1 Asymptotics of W^r^^^(A) 

In this subsection, we construct the modifiers J^j^^-^{X) and give the asymptotics of when 
A — > +00. For simplicity, we shall omit the lower index (-t-00) in all the objects defined hereafter. 

We first look at the problem at fixed energy (i.e. we take A = in the previous formulae). Hence 
we aim to construct modifiers with scalar phases (^^(x,^) and matrix- valued amplitudes p^(a;,^) 
such that the amplitudes c^(a;,^) of the operators = HJ^ — J^v^{Dx) be short-range in x when 
X +00. We adapt here to our case the treatment given by Gatel and Yafaev in [5] where a similar 
problem was considered in Minkowski spacetime (see also our recent paper [3|). 

The operators are clearly FIOs with phases ip^{x^S) and amplitudes 

c^{x,O^B^{x,i)p^{x,0-iT^d,p^{x,0, (3.19) 

where 

B^{x, = T^d^^^(x, e) + a(a;)r2 + &(x)r" + c(x) - (3.20) 

As usual we look for phases ip^ close to and amplitudes close to 1. So the term dxP^ in (|3.19p 
should be short-range et can be neglected in a first approximation. With = 1, we are thus led to solve 
B^ — 0. However a direct calculation leads then to matrix-valued phases (p^ whereas we look for scalar 
ones. We follow [3] and solve in fact [B^Y = 0. Using crucially the anticommutation properties of the 
Dirac matrices (|2.16p . we get the new equation 

[B^f = (d^ip^f +a^+b^ + ic- v^f + 2(c - v^)(B^ -c + v^)^Q. (3.21) 

If we put B^ = in (|3.2ip . we obtain the scalar equation 

r±(x, C) := {dxV^f +a^+b'~{c- v^f = 0. (3.22) 

We look for an approximate solution of (|3.22p of the form (/3^ (2;, = + 0^ (.t, ^) where 0^ (x, ^) should 
be a priori relatively small in the variable x. Recalling that (v^)'^ = + vr? by p. lip , we must then 
solve 

2^9x0^ + {dA'^f -f + {\? - m^) - + 2cv^ = 0. (3.23) 
If we neglect {dx(j)'^)'^ (|3.23p . we finally get 

2^dx<l}^ = -[a^ + {b^ - m^) - -f d^] , (3.24) 
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where we have introduced the notation d^{x,^) = 2c{x)v^{^). Note that by (|2.20p and (|3.1ip . the 
following estimate holds 

Va,^eN, \d^d^d^{x,0\ <Cap{x)-^-°'{£_y-'^, Vx e M+, G M*. (3.25) 

Therefore, using (12.201) again and the previous estimate ()3.25p . we see that — is short-range when 
X — > +00 whereas — ???? and are long-range (of Coulomb type) when x — s- +00. Hence we can define 
two solutions of (|3.24p for all ^ 7^ as follows 

'^^(^,0-2|X [^'(")-^'(^)]^"-2|i, [(^'W-"^')+^^(^'^)]^^+2|X m-m)'ds. (3.26) 

Remark 3.3 Let us emphasize that we only add the quantity ^ /o^°°(^(s) ~ m)^ds in i3.26\) in order to 
prove that the Isozaki-Kitada and the Dollard modifications coincide (see Theorem \3.!^) . In the general 
case however, the phases (f>^{x,S^), solutions of | j'. 24\ ) would clearly take the form for all ^ ^ 

(j)±^x,i) = — [a'{s)~c^{s)]ds- — {b'{s)-m^)ds ^ c{s)ds + C{^), (3.27) 



where C (^) is a constant of integration. 

With this choice, we obtain for ^ 7^ (see p. 221) ). 

r±(x,0 = {d.<t>^ f - ^ [a\x) + {b\x) - m2) - c\x) + d^{x,0]^ ■ (3.28) 

Moreover it is easy to see that the rests satisfy the estimates 

Va,/3eN, |a^afr±(x,OI <a/3(a;)-'-"(e)-^ Vx £ M+, G M*. (3.29) 

In our derivation of the phases (|3.26p . it is important to keep in mind that we didn't find an approx- 
imate solution of = but instead of (-B^)^ = 0. Therefore we cannot expect to take = 1 as a 
first approximation and we have to work a bit more. So we look for p^ such that B^p^ be as small as 
possible. According to (j3.2ip and ()3.22p . we first note that 

= r± +2(c- z.±)B±. (3.30) 

We find now a relation between B^ and {B^)^ . Using (|3.20p and (|3.24p . we can reexpress B^ as 

B^ ^ B^ +2iy^K^, (3.31) 

where 

= + (3.32) 

= j^{a^ + ib^ -m^)-c'' +d^)T^ +aT^ + ib~m,)T° + c]. (3.33) 

If we take the square of p.3ip we get 

{B^y = [Bff + 2v^B^K^ + 2v^K^B^. (3.34) 
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But from ((3?32| and (|3?TT|l we see that (B^)^ = -2v^B^. Whence (fS^Sil) becomes 

{B^f ^ -2v^B^{l - K^) + 2v^K^B^. (3.35) 
Now we replace the expression p.35p of [B^Y i^^to (|3.30p and we obtain 

r± = -2v^B^{l - K^) + 2v^{l + ~ (^■^^) 

We would like to isolate B^ in (|3.36p . We thus need to invert the functions {\ + K^ — -^). Using (|2.19|) . 
(|2.20p and (|3.25p . we get the following global asymptotics for 

Let us consider the set X = e R, |^| > R] where R » 1 is a constant. It follows immediately 
from the asymptotics p.37p and those of -^rj^ that (1 + — and (1 — K^) are invertible for all 
(x, ^) e R X X if the constant R is assumed to be large enough. In consequence we can write p.36p as 

i?±(l - = ^(1 + X± - ^)-V±(l - + (1 + X± - ^r'Bt. (3.38) 

for all {x,C) e R X X. The first term in the right hand side of ((3^ is small thanks to ([3^ but the 
second one is not. We choose in such a way that they cancel this term. To do this, we observe that 
the Fourier representations of the projections P™, i.e. the operators 

P±\i) = {tA-^^^'^ + "^r°)) = \ih± ;'f , (r^e + mr")) , ve^o, (3.39) 

+771 / /\ 1/4 +771^ ' 

satisiy the following equations 

i3o^(0^r(0 = 0, (3.40) 
by Lemma [5TT] and (|3.32p . According to (|3.38p . a natural choice for is thus 

P^ = {l-K^)-'P^{0. (3.41) 

for which we have 



q± := B V ^ 1 (1 + X± _ ^)-ir±(l - K^)-^P^{0. (3.42) 



Let us summarize the situation at this stage. For ^ 7^ 0, we have defined the phases ip^{x^^) = 
+ 4>^{x,£,) by l|3.26p and for ^ e X, the amplitudes are given by p.4ip . Directly from the 

definitions and from the asymptotics (|2.19p and (|2.20p of the potentials a, b, c, the following estimates 

hold. 

Lemma 3.2 (Estimates on the phases, the ampUtudes and related quantities) For all x G R+ 

and ^ G X with R large enough, we have 

V/3 e N, |af(/.±(x,OI < C;3log(x) (e)-^. (3.43) 
V|a| > l,V/3 e N, \d:dU^{x,i)\ < C^pixy^iO'^. (3.44) 
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\dlM(.x,0-xO\<§^- (3.45) 

Va,/3eN, |9^5fi^±(x,e)| <C„0(x)-i-"(e)-^-^. (3.46) 

Va,/3eN, \dSd^^{p^ix,0-P^{0)\<C^p{xr'-''{0-'-''- (3.47) 

Va,/3eN, |9^9fr±(:r,C)<C„0(x)-2-"(^)-^. (3.48) 

Va,/?eN, |5^5fg±(a;,e) <CaM^>"'""(0''"''- (3-49) 

Va,/5 e N, |9^9fc±(x,0 < C^a/3(x)-2-"(0-i-'5. (3.50) 

Thanks to p.43p . (I3.44|) . p.45p and p.47|) . for R large enough, we can define precisely our modifiers 
as bounded operators on H (see [25j for instance). Let G C°°(K) be a cutoff function in space 
variables such that x^(x) = if x < 5 and x'^i^) = 1 if a; > 1. Let also 9 G C°°(R) be a cutoff function 
in energy variables such that 9{£_) if |^| < ^ and 6'(^) = 1 if |^| > 1. For R large enough, are the 
Fourier Integral Operators with phases ip^{x,£^) and amplitudes 

P^{x,O=X^ix)p^ix,O0{^). (3.51) 

We finish this part by a first application of the previous construction. In the next Theorem, the 
modifiers J^^^ are shown to be time-independent modifications of Isozaki-Kitada type equivalent to the 
DoUard modification (|2.21|) . Precisely we have 

Theorem 3.3 For any ijj €z Ti such that supp G X , we have 

- ^ hm^e"«J± ^)e-*''"(^^)p£^V. (3.52) 
Proof: We only sketch the proof for the case (+). By definition of P™, we have 

'"/.'[("" "(" )]ds 

[/(t)P™V' = e"''*^^^"^e V°i+"' " ■ := V"(t)P+ V'- (3.53) 

Then, we write : 

e'*^J(+^)e-"''^(^-^P^V = e^*^y(t)(F*(t)e-^*''''(^="))(e'*'^^(^-'j+^)e-^*''"'(^-))P_7'V (3.54) 

= e'*"V{t) e'/ol-l-i- (e**''^(-°=")j+_^)e-^*''^(-°-))P^V'. (3.55) 
The classical flow associated with the Hamiltonian v^iC) — sgn{S)\/^^^^~\^Tn? is given by 

^\x,0 = {^ + t^==,0- (3.56) 

Then, using Egorov's theorem, we see that (e**"^'^^"'' Jj^^^e"**"^^^"^) is a FIO with phase ip+{t,x,^) = 
x^ + (/)+(x + trj, ^), and with principal symboQ P^{x + tr], ^) where 77 = — ^ ^ ^ 

+ TO^ 



^It means that the others terms of the symbol are o(l) when t — > +co. 
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Thus, e'/ol-l'*'' (e'*''^('°-)j+_^^e-**''^ is a FIO with the same principal symbol and with phase 
iPi(t,x,S,) = + (j)i{t,x,^) where 



-1 r + oo -1 px+tr) 

^tit.x.O^TT. / [a^(.)-c2(.)Ks-- / [{b'{s)-m-') + 2c{s)u+{0)]ds 

Jx+tri Jo 



tr, 



+ ^ (Ks) - mfds + [(&(sr;) - + ci.sr,)]ds (3.57) 

Since — / — c^(s)]c?s = o(l) when t +00, and by making a change of variables in the last 

Jx+t-n 



integral, we obtain 



1 rx+tri 1 

0^(i,^,O=-^y^ [(62(s)-™2) + 2c(s);.+ (0)]ds+- / {b{s)-mfds 



1 



2(6(s) - m)m + 2c(s)i^+(0]ds + oil). (3.58) 



Using again that / [(fe^(s) — m^) + 2c(s)i^+(^))] ds = o(l), we see that 

Jtrj 

1 /•*'? 1 /•+00 

^tit.^^O^-Yj^ [{b^{s)-m^) + 2c{sy+{0)]ds + -J^ {b{s)-mfds 



1 



t»7 

2(6(s) - m)m + 2c(s)i^+(0]ds + o(l). (3.59) 



Then, 

cj,+ {t,x,0 = -7r, iKs)-mfds + +- {b{s)-mfds + o{l)^o{l). (3.60) 
Jo Jo 

Using (I031) . ([04]) . (IOT]) and the continuity of FIOs, we see that 

^^J^[...]ds l^e'^'^iD^) j+_^^e-^t^^{D.)-jp^^ = p^^ + o{l) (3.61) 
and Theorem 3.3 follows from (|3.55p and p.6ip . (} 

We now construct the modifiers at high energy J^^-fW so that they satisfy p.9p and p.lSp . We 
still omit the lower index (+00) in the next notations. Comparing p.lSp and p.52p suggests to construct 
J^(A) close to e"*"**^ J='=e*'^^ which are clearly FIOs with phases ip^{x,£^,X) = + 0^(a;,^ + A) and 
amplitudes P'^{x,£^ + A). 

With J±(A) = e"'^^J±e*^^, we see from ([XSUI that the amplitudes 

c± (x, A) - B± (.T, C + A)P± (.X, C + A) - zri9,P± (x, e + A), 

of the operators C*(A) — H{\)J^ {X) — J^{X)iy^{Dx + A) would satisfy the estimate 

c±(x,e,A) =0((a;)-2^"i), (3.62) 

for ^ in a compact set. Here and in the following, the notation f{x, A) = 0{{x)~'^X~^) means that f{x, A) 
decays as (a;)~^ when x +00 and as A""'^ when A — > +00. We want however the amplitudes c^{x, ^, A) 
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to be of order 0((a;)~^A~^) and the decay in (|3.62p is not sufficient for our purpose. In consequence, we 
need to refine our construction. Following the procedure given in [4 , we look for modifiers J^(A) defined 
as FIOs with phases ^, A) and with new amplitudes P^{x, ^, A) that take the form 



P^{x, e, A) = [p^{x, e + A) + e + X)l^{x) + ^P^k^ix) 



(3.63) 



(up to suitable cutoff functions defined later), where P± denote the projections onto the positive and 
negative spectrum of F^. Here the correctors l^,k^ (that can be matrix- valued) will be functions of x 
only and should satisfy some decay in x (see below). It will be clear in the next calculations why we add 
such correctors to the amplitudes p^{x, C + A). 

We now choose and in p.63p so that the amplitudes 



c±(x,e,A) 



i?± (a:, e + A) [p± ix,^ + X) + jp^ix,^ + X)l^{x) + ^P^k^ (x) 



(3.64) 



d^p^ (x, C + A) + jd^p^ {x,^ + X)l^ (x) + ip± (:e, ^ + X)dJ^ (x) + ^P^d^k^ (x) 



of the operators C*(A) be of order 0((a;)^^A^^). 

To prove this, we need the asymptotics of the different functions appearing in (|3.64p . For x in 
and for A large enough, we obtain (after long and tedious calculations) 



r.±(e + A) = ±[A + e+^]+ OiX-% 

d±(x,e + A)=±2c(x)[A + e+^] + 0{{x)-^X-'). 

K^ix, e + A) = ±^[2P^ c{x) + a{x)T^ + (&(x) - m)r°] + 0((x)-i A""). 

P±{^ + X) = P±+ O(A-i). 
P^{x,^ + X)=P±+ O(A-i). 

9.p±(x,e + A) = ±^P^(a'(x)F2 + 6'(x)F0) + 0{{x)-'X"^). 

B^ix, ^ + A) = T2(e + X)P^ + 2c{x)P^ + a{x)T^ + h{x)V^ + 0{X-^). 

g±(a;,^ + A) =i?±(x,e + A)p±(x,e + A) =±i^c2(x)P±+ 0{{x)-^X-^). 



(3.65) 

(3.66) 

(3.67) 

(3.68) 
(3.69) 

(3.70) 

(3.71) 

(3.72) 



We mention that the following simple equalities have been used several times to prove the preceding 
asymptotics 



1 + ri 



h 




2P+ 



1-Fi 




h 



2P_ 



(3.73) 



By ((3?69ll . ((370)) . |3JT|l and ((372)) . the amplitudes c=^(a;,^, A) take the form 
c±(x,C,A) = ±T^c2p± i^^. 



±^c^P± ± -^c^P±l^ 



1 

■2A 
1 

"A2 
-zFi 
1 

"a 



T2(e + A)P:p + 2cP^ +aF2 



6F" 



0{ 



P^k^ 



± ^P^{a'T^ + 5'rO) ± ^P^[a'T^ + b'Vy^ 



P±+ 0{j))dj^ 



^p^d..k± 



1 



(.t)2A2 
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From the asymptotics (|2.20p of the potentials a, 6, c, we rewrite this last expression as 

c±(a;,^,A) =±^c2p± T ^^t^* T ^r^PTi^T^ + bT'') - jT^P±dJ^ + R{x,X), (3.74) 
where the rest R{x, A) satisfies 

^(-' ^) - Q( (,)2;,2 + + + + j • (3-75) 

Now we choose the correctors l^, in such a way that the terms of orders 0(A^^) in (|3.74[) cancel. Once 
it is done we shall have to check that the rest (I3.75P be of order 0((a;)~^ A^^). 

There are clearly two different types of terms in the expression (j3.74p : on one hand the terms 

2A A A^2 ^' 

"live" in Ti.± = P±{'H); on the other hand the terms 

TjP^k^ T^rip^(a'r2 + 6'rO) = ip^[T2fc± + '-{a'T^ + bT')], 

"live" in TYip = Pzp{'H). Since the Hilbert spaces 7i_ and H+ form a direct sum of H, i.e. 7i = H- ® 
we can consider separatly the equations 

±^c^TidJ^ = 0, (3.76) 

T 2fc± + ^(a'r^ + 6'r") = 0, (3.77) 
in order to cancel the terms of order 0(A^^) in (I3.74p . We solve first (|3.76p and obtain 

• r + oo 

l^{x) = l{x) = - / c^{s)ds. (3.78) 

^ Jx 

Then we solve (I3.77P and get 

fc±(x) = ±^(a'(x)r2 + b\x)T'>). (3.79) 
The functions ^ and clearly satisfy when x — > +00 

l{x) ^ 0{{x)-^), dj{x)=0{{x)-^), k^{x)^0{{x)-^). (3.80) 
Finally with this choice of correcting terms I and k^ , we conclude from (|3.74p and (|3.75p that 

c=^(a;,C,A) = R{x,X) = 0{{x)-^\-^). 
In fact, we can prove that for all x E M+, ^ in a compact set and A large enough 

Va, f3eN, \dy^c^{x, A)| < C^^ (x)-^-" A-^. (3.8I) 

Let us summarize the previous results. The modifiers J^(A) are (formally) constructed as FIOs with 
phases (p^ {x, ^, A) = + 0''^ (x, ^ + A) where 

0±(a:,e + A) = ^^^^ (^y [a\s)~c\s)]d.s- [{b^i-s) - m') + d^i-s,^ + X) ds] 



21 



+ J {b{s) ~ mydsj , (3.82) 

and amplitudes 

P±(x,^,A)= [p±(a;,e + A) + ip±(.T,C + A)Z(a;) + ^P^fc±(x)], (3.83) 
where I and are given by (|3.78|) and (|3.79p respectively. 

Unfortunately, since cj)^{x,^ + A) = 0{< x >) when x — oo, this phase does not belong to a good 
class of oscillating symbols. So, we have to introduce some technical cutoff functions in the amplitude 
in order to localize x far away from —oo. Moreover, these cutoff functions must be negligible in the 
asymptotics in the previous calculus. We follow the strategy exposed in [22] which we briefly recall here. 

We consider a fixed test function ip G C^(M) and we want to calculate the asymptotics of VV^^^{X)Tp. 

Since -0 ^ C^(M), at high energies, translation of wave packets does not dominate over spreading. So we 
introduce a cutoff function (depending on A) in order to control the spreading. 

Let xo G C'^(IR) be a cutoff function such that Xo(0 = 1 if U l< 1> Xo(0 if | ^ |> 2. Using the 
Fourier representation, we have easily : 

Ve > 0,V7V > 1 , II [xo(^) - 1] IIl^(e)- 0(A-^). (3.84) 
Now, let us define the classical propagation zone : 

n = {x + t ; xe supp tp, t e K+l, (3.85) 

and let 77+ S C°°(R) be a cutoff function such that 77+ = 1 in a neighborhood of fl and 77+ = in a 
neighborhood of —00. We consider 

^^(A) = (77+ - 1) e-"''*(^^+^) P^^^^,(^£^)^ (3.86) 

Liemma 3.3 For A > 1, e e]0, 1[, t e M.^, and N>1, we have : 

||X±(A) |U.(R)=0(<t>-^A-^). (3.87) 

Proof: We only sketch the proof for the case (+). Using the Fourier transform and (|3.39p . we easily see 
that 

K-iX) ^ I i,Hx) J(J e-^^) (a + ^)||=|^) XoiO di) m dy. (3.88) 



where (^(0 = A'(x - y)£, - t^J {\^^ + Xf + in"^ . So, 

Since ^ is in a compact set, e < 1, y G supp ip, we easily obtain for x G supp (77+ — 1), and X >> 1, 

\d^^{0\>cX'{l + t), (3.90) 
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for a suitable constant c > 0. We conclude by a standard non stationary phase argument. <^ 

Now, we can define precisely ours modifiers J* (A) in order to calculate the asymptotics of Wi^ ^ (A)-;/'. 



According to (|3.84p . it suffices to calculate the asymptotics of W^^~^{X)xo{^^)'4'- We first remark that 

for A ^ f and e < f , we have S, + X£Xif-^£ supp xo- So, we can define the modifiers J^(A) as FIOs 

A^ 

with phases 1^3^(2;, A) = + (f)^{x, ^ + A) where 0^(a;, C + given by p.82p and with amplitudes 



P±(x,^,A)=77+(x) [p±(a:,e + A) + ip±(x,e + A);(x) + ^P^A;±(x)] xo(^), 

where I and are given by p.78p and (|3.79p respectively. 

With this definition, we can mimick the proof of Theorem 13. 31 to get 
Lemma 3.4 For ifj e C^(R) and for A large, we have 



(3.9f) 



W, 



,D 



(+00) 



(3.92) 



Moreover, it is easy to see that the estimates p. 811) are still satisfied, so we can prove our main 
estimate p.9p . Precisely we get 

Lemma 3.5 For ijj G C|^(R) and when A tends to infinity, the following estimate holds: 

Proof: Everything done in [4l Lemma 3.3 works here in the same way. All the contributions coming 
from the cut-off function 77"*" are negligible using the same arguments as in Lemma 13.31 since the support 
of the derivatives of rj^ are far away from f2. <^ 



We end up this section giving the asymptotics of W^(+oo) when A is large. According to Lemma[i 
we have for any ip € C^(M; C), W^(+^)(A)V' = 0{X-^). Thus we only need to compute the 

asymptotics of the modifier jfl^^oo) (•^) shall consider as pseudodifferential operators with symbols 

j±(x,C,A)=e***(-'«+^)p±(x,^,A). 
Using the explicit expressions (|3.82p and (|3.9ip . we first get the asymptotics 



'-±(x,^ + A) =T / c{s)ds + j^ 



+ 00 



+ / {b{s)~myds + 0{^^) 



(a^ - c'){s)ds - I {b^{s)-m^)ds 
log(x) 



A2 



P^{x,^,X)=v+ix) 



P± ± ^P^ (aF^ + fepO) 



A 



Moreover using a Taylor expansion of e* at t = 0, we get from (|3.93p 



(3.93) 
(3.94) 

(3.95) 
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with 



C+{x)= c{s)ds, C+[x)^ [a^ ~ c^){s)ds - [b'' [s) - m^)ds + [b{s) - mfds. {iM) 



Combining now (|3.94p and p.95p . we obtain 



i A) = e 



l{x) 



P± +^C+(x)P± ±-P^{aT^+hT^)+ ^P± + O(-). (3.97) 



2A 



+ CXD 



But notice from (P?75)) that 

M= _(^y^ a\s)ds- {b%s) - m')ds + {b{s) - mfds 

and from the anticommutation properties (|2.16p of the Dirac matrices that 

P^iaT^ + bT°) = {aT^ + bT")P±. 

Hence p.97p becomes 



j±(:E,C,A) =eT*^^(")77+(x) l + 7^f/ a\s)ds- {b^{s) - m^)ds + {b{s) ~ mfds 



Eventually, if we introduce the notations 



= a^{s)ds~ f {b^{s)~m^)ds+ f (6(s) - m)^^^) ± 1 (aF^ + 6r°), 



'^^ Jx Jo Jo 

we deduce from ((XM)) and the fact that 77+ (x) = 1 on supp ip, the following Proposition 
Proposition 3.1 For any tJj e Cf^iR;^^), 

where {x) and [x) are given by S3.96\) and S3. 99\) respectively. 



(3.98) 



(3.99) 



(3.100) 



3.2 Asymptotics of W(too)('^) 

In this subsection, we focus on what happens at the event horizon and give the asymptotics of W^^s^ (A) 
when A +00. In fact, we shall derive them from the results obtained in the preceding subsection 13.1 
after some simplifications of our model. As usual, we shall omit the lower index (—00) in the objects 
defined or used hereafter. 

Recall that the expressions of the wave operators at the event horizon are given by (see (|2.22p ) 



W^ = s- lim e'*"e-'*"°P^, 

t->±oo + 
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where Hq = T^D^ + cq, H = T^Dx + aT^ + toF^ + c and the potentials a,b,c — cq satisfy (|2.19p when 
X — > — oo. We first simphfy this expression in a convenient way. Let us introduce the unitary transform 
U onH 

jj ^^-^^'c-(x)^ ^-(x)^/ [c{s)-CQ]ds+ cox, (3.101) 



and define the selfadjoint operators on Ti 

Ao^r^Dx, A^U*HU. (3.102) 
Using (|3.10ip . a short calculation shows that the operator A can be rewritten as 

A = r^Dx + W{x), (3.103) 

where 

Wix) = e'^"^'<'''\aix)T^ + bix)T°)e-'^"^'^''\ (3.104) 

Note that according to the anticommutation properties (|2.16p of the Dirac matrices, the potential W 
satisfies WT^ + T^W = and W'^{x) = a'^{x) + b'^{x). Moreover from ((2?T9l) . we get the following 
estimates for W 

3a >0, M^(a;) = 0(e""^), a: ^ -oo. (3.105) 
Using the unitarity of U and (|3.102p we rewrite as 

= Us- lim e^^'^U*e~'*"°P^, 

= Us- lim e'^A^-^tA„^^tAo^*^-itHop (3.106) 

Now we can simplify the strong limit appearing in (|3.106p in two steps. First we claim that 

s-^lim e'*-^''[/*e-**-f^"P:p =e*r'"°^P:p. (3.107) 
Indeed, using the particular diagonal form of F^ given in (|2.18p and since e~'*-f^o = g-^t^og-jtco^ have 

^itAajj*^-itHop_^ _ gitAo^ir'-C^(x)^-iAo^-itcop_^ _ ^ir'- (x^t) ^-itco ^3 ^08) 

When t +00, the right-hand-side of p.lOSp can be written using (|3.10ip as 

^-,C-ix-t)^-^tcop_ ^ ^-^{f-^'('^('^)-'^o)ds+cox^ 

from which (|3.107p follows when t — > -t-oo. The case t — > — cx) is obtained similarly. 

Second since the potential W decays exponentially when x — oo by (|3.105p . it follows from the 
methods used in [31 [H] that the wave operators 

W^{A,Ao) = s- lim e**^e-**^°P:n=, (3.109) 

t— >±oo 

exist on TC. Hence by p.l06p . (|3.107p . (|3.109p and the chain-rule, we obtain the following nice expressions 
for 

= UW'^{A,Ao)e'^''=°'' P^. (3.110) 
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At last since U and e*'" '^"^ commute with e'^^, it is clear from (|3.110|) that it is enough to know the 
asymptotics of 

W^iA,Ao,\) ^ e-'^'''W^{A,Ao)e'^'-' 

when A —> +00 in order to get the asymptotics of W^{X). 

Note here that the A-shifted wave operator W^{A, Aq, A) is exactly the kind of wave operator studied 
in our previous paper [4] in which the asymptotics of W^{A, Aq, A) were calculated. Nevertheless we can 
also easily derive these asymptotics from the results of the preceding section. For completeness this is 
what we choose to do here. 

We thus follow our usual strategy and construct modifiers J(f(A) corresponding to W^{A, Aq, X). 
This problem is in fact similar to the one in subsection 13. II It suffices to replace iJ™ by Aq and H by 
A in our calculations. From the explicit form (|3.102p and (|3.103p of the operators Aq and A, we deduce 
that we can use the results obtained in subsection 13.11 with the following changes: (1) Since the mass 
m doesn't appear in Aq hence we take m = 0. (2) The long-range matrix-valued potential b and scalar 
potential c don't appear in A (see p.lOSp and (|3.105|) ) hence we put b{x) = c{x) = 0. (3) The short-range 
matrix- valued potential a{x)r'^ is replaced by W{x). (4) The projections P™ are replaced by Pzp since 
we work at the event horizon. Noting that these changes also entail that — ^nd d^{x,£^) — 0, 

we obtain the following results. 

At fixed energy A = 0, the modifiers are defined as FIOs with phases 



and amplitude^ 



1 r W^(t'\ 1 

p^{x,0 = {l~K^{x,Or'PT^ K^ix:0=T7r. - -^'^' + W{x) . (3.111) 



2C 



2^ 



ip'^{x,^,X)=x^+—^-— I W'^{s)ds, (3.112) 



At high energy, the modifiers J^{X) are defined as FIOs with phases 

1 

2(?+A) 
and amplitudes 

P^ix, A) = p^{x, e + A) + ^P±k^{x), (3.113) 

where k^{x) = ^jW'{x). Using these definitions and p.lOSp . we can prove that the symbols c^(a;,^, A) 
of the operators C^(A) = A{X)Jq{X) — J'^{X)Aq{X) satisfy the estimates 

Va,/3eN, |9^afc±(x,C,A)| <C„^ (3.114) 

for all x £ and A large enough. Finally as in the proof of Lemma [3.51 the estimates (I3.114p are the 
main ingredients to prove the equivalent properties to (I3.14p and p.9p . Precisely we have 

Lemma 3.6 For any ip G C^(M;C*) and for X large, the following estimate holds 

\\{W^{A,Ao,X) - = OiX-^). 



^In the same way as the preceding section, we should add some technical cutofT functions which are negligible in the 
asymptotics. 
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We now use Lemma [3.61 to compute the asymptotics of VF^(yl, Aq, A)^ up to the order 0(A ^). For 



any e 



and for A large, we have 



W^iA,AQ,X)ij^4{X)ij+ O(^). 

Hence, it is enough to compute the asymptotics of JoiX) for A large. Using p.llip . (|3.112p . p.ll3p and 
after some calculations, we obtain 



Jo±(A)^ 



1 + 



2A' 



(3.115) 



Note that we retrieve naturally the same formulae as in [4]. Eventually combining (I3.110p and (13.115^ . 
we obtain the asymptotics of W^{X) for A large 

Proposition 3.2 For any -0 e C^(K), 

W^_^^{X)i:^u[l + jQ^{x)\e^^'-"'-'P^^+ O(^), (3.116) 

where U is given by IIS.101\) . Q±(x) = i J^°° W^{s)ds T W{x)^ and W{x) is given by 19.1 04\ ). 

3.3 Proofs of Theorems [SH] and [372] 

In this last subsection, we use the asymptotics of W^^^{X) obtained in Propositions 13 . 1 1 and 13 . 21 to prove 
the reconstruction formulae given in Theorem 13.21 and finally prove Theorem [371] 



Proof of Theorem \3.2\ We only treat the case of the transmission operator Tr and give the proof of (|3.ip 
since the proof of ()3.2p corresponding to the transmission operator Tl is similar. Recall that we want to 
compute the asymptotic expansion when A +oo of 

Fi{X) =< rfle'^^V,e^^> >=< M/(+o,)(A)V, W^(+oo)(A)0 >, 

for & C^(M;C'^). Using Propositions 13.11 and 13.21 and the notations therein, we have 



FiiX) = <e 



iC+ix) 



P-iIj,U 



l + ^Q+{x) 



iV^ cqx 



P-cb > + O(^), 



> 



> 



1 



+ 0(^). 

We now compute separatly the terms of different orders in p.ll7p . 
Order 0: Since r^P_ — — P_, the term of order reads 



(3.117) 



(3.118) 
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Moreover from (j3.96p and (|3.10ip , the phase C+ (x) — C (x) + cq takes the simple form 

.0 

C+{x) - C^ix) + cox = - / [c(s) - co]ds + cqx. (3.119) 

<y —oo 

Order 1: Using r^P_ = — P_ again, the term of order 1 can be written as 

^ ^t[C+{x)^C-{x)+coa:] _^ {Q+ )* ) P_ ^/j , P_ (f) > . 

Since W'^ = a"^ + 6^ and WP- = e'^'^' {aV'^ + 6r°)P_ by (|2l6| . the term {Q+)*P- takes the form 

iQ^yP- = {-\j (a'+6')(s)ds--e2'^"(ar2 + 6rO)jP_. (3.120) 
Moreover from (|3.99p the term i?^ is 

i?- = M a^{s)ds-- ib^is)-m^)ds + - ib{s) - m)^ds - -{aT^ + bT°). (3.121) 

2 Ja; 2 7q 2 Jq 2 

Hence adding (|3.120p and p.l2ip . the term of order 1 reads 

Ac-^ix)-c-ix)+cox]n a^(^s)ds + ^ f b\s)ds+'- {b{s) ~ mfds + ^m^x)p^ij,P_cP > 

»[C+(:r)-C-(:.)+Cox]|^lg2»C-(^p2 ^^pO) _^ 1 (af^ + bP")) P_ ^/;, P_ ?i > . (3.122) 

FinaUy using that e'[<^+(^)-c^"(^)+'=o=^l is scalar, that {aV^ + bT°)P± = P^{aT^ + bV") by (^1^ and the 
fact that < P+ijjjP-'t' >— 0; we see that the last term in p.l22l) cancel, i.e. 



< e 



< e 



[c+(x)-c-(x)+co.] (ie^'C- (ar2 + 6r") + ^(ar^ + 6r"))F_^, P_(/. 0. 



Hence the term of order 1 is 

- a\s)ds + - h\s)ds + - {b{s)^mfds + -n?x)P-^,P-^>. 

2 J-oo 2 2 Jq 2 / 

(3.123) 

If we introduce the following functions 

0(x) — g-*/-=„[c(s)-co]<is+ico3: 

/• + 00 ^.0 /•+00 

A{x) = e(x)( / a2(s)ds+ / b^{s)ds+ / (fo(s)-m)2(is + 771^2;), 

^J-oo J -co Jo ' 

we have proved the reconstruction formula p.ip and thus Theorem 13.21 ^ 

Proof of Theorem \3.1\ We show here that the reconstruction formula (j3.ip entails the uniqueness of the 
parameters M and Q under the additional assumption that the charge q of Dirac fields is known, fixed 
and nonzero. The same result can be shown from the reconstruction formula (|3.2p in a similar way. 
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We first compute one of the integrals that appear in (|3.ip which will be useful in the later analysis. 
Using the explicit expressions of F, ai given in (|2.2|) and p.isp as well as the definition of the Regge- 
Wheeler variable x{r) given in (|2.6p . an easy calculation shows that 

af{s)ds^il + h'-, (3.124) 
2 ro 

where ro is the radius of the event horizon. 

Now let us consider two transmission operators T/^i and corresponding respectively to parameters 
Mj,Qj,mj,{j = 1,2) and qi = q2 = q where q is supposed to be known and nonzero. In what follows, 
all the objects corresponding to T/.j with j — 1,2 will be denoted by the usual notations with a lower 
index j. We suppose that T/^ — Ti^2- In consequence we also have i^;,i(A) = F/^2(-^)- Our goal is to 
prove that Mi = M2 and Qi ~ Q2- Using Theorem 13.21 and identifying the terms of same orders in the 
reconstruction formula ()3.ip . we thus get 



ei(a;) = 62(0;), (3.125) 
Ai{x) = A2{x). (3.126) 

By p.3p and a standard continuity argument, ()3.125p leads to the equality 

/O M 
[ci(s) - co,i]ds + ico^a; = / [c2{s) - co^2]ds + ica,2X + 2fc7r, (3.127) 

where fc G Z. If we derivate p.l27p with respect to x, we obtain 

co,i = co,2 := Co. (3.128) 
Now by p.l24p . (|3.126p leads to the equality 

{l + -f + - / bl{s)ds + - / (61 (s) - mfds + -mix 

2 ro,i 2 2 Jo 2 

= {1 + 7^f— + 7; / bl{s)ds + - / (62(5) - mfds + -mix (3.129) 
2 7-0,2 2 J^^ 2 Jo 2 

If we derivate (|3.129p with respect to x, we first get 

mi = TO2 := 7n. (3.130) 

Hence the mass m of Dirac fields is uniquely determined. Moreover, using p.l30p . (|3.124p and the 
homogeneity in the parameter I, we obtain from (j3.129p 

^0,1 = ''0,2 := ro- (3.131) 

Therefore the radius ro of the event horizon is also uniquely determined. Now if we combine p.l3ip and 
Co = ^ into (|3.128p . we get (since q is supposed to be nonzero) 

The charge Q of the black hole is thus uniquely determined. Eventually since ro cancels the function F ^ 
we get from (|2.2p that 

Ml = M2 :=A^- ^ii^, 
2ro ' 

and the mass AI of the black hole is uniquely determined. This finishes the proof of Theorem 13.11 <^ 
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4 The inverse problem for dS-RN black holes (A > 0) 



In this Section, we study the inverse problem in the case A > corresponding to dS-RN black holes. 
In a first part, we prove the same kind of results as in Section [31 that is we prove that the parameters 
M, Q and A are uniquely determined by the high energies of the transmission operators or Tr. In a 
second part, we prove by means of a purely stationary method that the parameters M, Q and A can also 
be uniquely determined from the knowledge of the reflection operators L or R on any interval of energy. 

4.1 The inverse problem at high energy 

As in Section [31 we shall assume here that one of the following functions of A G M 

Fi{X) =< rfle^^^V,e^^> >, G/(A) rie^^^V,e^^> >, 

is known for all large values of A, for alH e N and for all ^, e with V': G C'5^(R; C^). We emphasize 
that in this case the construction of the modifiers are simpler than in the previous section due to the 
decay of the potentials at infinity; the phases of the modifiers constructed later will belong to a good 
class of oscillating symbols. In particular, we do not need a technical cutoff function ry+ and a cutoff 
function xo in order to control the spreading of the wave packets as in Section [3] and we can consider test 
functions & Ti- with tjj,(j> G C(f (R;C^). We also assume that the mass m and the charge q of the 
Dirac fields are known and fixed. Furthermore the charge q is supposed to be nonzero. Then our main 
result is 

Theorem 4.1 Under the previous assumptions, the parameters M, Q and A of the dS-RN black hole are 
uniquely determined. 

This Theorem will follow from the following reconstruction formulae obtained on each spin-weighted 
spherical harmonics 

Theorem 4.2 (Reconstruction formulae) Let tlj,(l) E H such that tp,(l) E C^(M;C'*). Then for A 
large, we have 

Fi{\) = <e{x)P-tl:,P^(l)>+\<A{x)P_i:,P_(l)>+0{X~^), (4.1) 

A 

Gi{X) = <e(x)P+7/',P+0> -i <^(x)P+^A,P+0> + 0(A-2), (4.2) 
where 6{x) and A{x) are multiplication operators given by 

e(x) = e-^^-^(=+-=°)^ A{x)^-[J {af{s) + b\s))ds)Q{x), (4.3) 
and a constant (3 given by 

(c(s)-co)ds+ / {c{s)-c+)ds. 

'OO Jo 

We shall prove Theorem 14.21 using the same global strategy as in the proof of Theorem 13.21 From (|2.30p . 
(|2.31[) . (|2.33p and the fact that e'^^^ corresponds to a translation by A in momentum space, we express 
F{X) and G{X) as follows 

Fi{X) = <M/-^)(A)^,I^+^)(A)0>, (4.4) 
Gi{X) = <WrJX)i:,W+iX)cj,>, (4.5) 
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with 

W± ,(A) = e-'^'^Wt ,e*^^= s- lim ^''^W e-'^"°'^^^ P^, (4.6) 
VF± - e^^^^W^i ,e^^^== s- lim e**-f^(-^)e-"-^+(^)p±, (4.7) 

( + oo)V ' (+00) t^±CX3 ^ ^ 

and 

H{\) = T\D, + A) + a{x)T^ + &(s)r° + c(x), 
Ho{X) = T\D, + A) + CO, H+{\) = T\D, + A) + c+. 
In consequence, it is enough to obtain an asymptotic expansion of the A-shifted wave operators W^^^s^ (A) 



in order to prove the reconstruction formulae (j4.1|) and (|4.2 

Note first that the A-shifted wave operators W^^.^{X) given by (|4.6p are exactly the same as in the 
case A = studied in Subsection 13.21 For completeness we recall here the asymptotic expansion of 
W^_^^{X) obtained in Proposition [3?1 For any ■tp € H, ■ip € Cg°{R; C*), we have 



W^-^^iX)^ = u[l + jQ^ix)y^"^'>-P^4'+ O(^), (4.8) 



where 

U = e-'^'^'^''\ C-{x)^ I [c(s)-co]ds+ Cox, (4.9) 



1 /• — oo 

Q±(x) = ^[ij W^is)ds T W{x)j, W{x) = e''^''^"(")(a(x)r2 +6(a;)rO)e-'^'^"("\ (4.10) 

Note second that the A-shifted wave operators W^^^{X) given by (|4.7p are very similar to (|4.6p . the 
constant cq being replaced by c+ and the projections P:^ being replaced by P± since we work now at the 
cosmological horizon. Hence they can be studied exactly the same way as in Section 13.21 Since there 
are slight modifications in some formulae, we recall here the procedure but omit the proofs. Using the 
unitary transform (|4.9p . we simplify the wave operators W^^^ as follows 

M/± ,=Us- lim e^*^e~"^"e**^''C/*e^*^+P±, (4.11) 

where we have used again the notations Aq = T^Dx and A = U*HU = T^^D^ + W{x) from (|3.102p and 
p.l03p with the potential W given by (|4.10p . We also recall that by (|2.16p this new potential W{x) 
satisfies the properties 

T^W + WT^ = Q, W'^^a^ + b^. (4.12) 

as well as the global estimate 

3a >0, W^(x) = 0(e-"l^l), Vx G M. (4.13) 

The potential W is thus very short-range both at the event horizon and at the cosmological horizon. Now 
an easy calculation shows that (to be compared with (|3.107p and its proof) 

s- lim e**-^«C/*e-*^+P± = e*^'^e*^''=+^P±, (4.14) 

>±oo 

where the constant f3 is given by 

(3=1 {c{s)~co)ds+ {c{s)-c+)ds. (4.15) 
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Furthermore it is immediate from (|4.13p that the wave operators ^o) = s— lim e'^'^e exist 

t — >±oo 

on Ti.. Hence we conclude by the chain-rule that W^^^ take the nice form (to be compared to the 



expressions (ISHOll obtained for W^_^^) 



W, 



(+00) 



(4.16) 



Since U and e''"^^e''"^^+^ commute with e'^^, we finally get the following expression for 



where 



Clearly it is enough to know the asymptotics of W^{A^Aq,X)P± when A ~\-oo in order to get the 
asymptotics of W^^-^[X). In fact the calculations are exactly the same to what has been done in 
subsection 13.21 (it suffices to replace by P± in these calculations) or in [3]. Hence we only give the 
final result without more details. For any -0 € Ti, -0 € C5"(IR; C''), we finally obtain 



1 ^- 



P±^+ O(^) 



(4.17) 



Proof of Theorem \4.^ We now use the asymptotic expansions (14. 8p and (I4.17P to prove the reconstruction 
formulae (|4.ip and (|4.2p . Since the proofs are analogous, we only treat (|4.ip . Using the previous notations 
we clearly have 



where U is given by gH), Q^{x) ^ ^(i /+°° W^{s)ds ± W{x)^ and W is given by glHl)- 



F,(A) =<u\l + |g-(x)le*i^'Vr'^+"F_V, [1 + ^Q+{x)\e'^''''"'P-^ > + 0(-^' 
L A J L A J A-^ 

Since U is unitary and since r^P_ = —P-, we reexpress (|4.18p as 

F{iX) = < e~'''~'("+~'=°)^P_7/;,P_0 > 

+ j < e-"3-'(^+-^°)"(g-(:r) + (g+)*(a:))F_0, > + O(^). 

From the explicit expressions of (5+ and , (|4.19p becomes 



1 



1 



(4.18) 



(4.19) 



(4.20) 



1 /J 
A^ V2 



(-y w^2(^)^^_^(^)jp_^^^p_^>^o(^). 



Eventually observe that W{x)P- — P+W{x) by (|2.16p and that < P^^,P-(j) >— 0. Hence we obtain for 

Km 



Fi{X) = < e-''^-'(^+-^")^P_0,P_0 > 

+ 7TV / VK2(s)ds < e-"3-'(=+-=°)"P_0,P_0 > + 0( — ). 
2A 



(4.21) 



32 



Denoting 

• p + oo ■ p+oo 

e(:E)=e-'5-(-+--»)-, A{x)^-[J W\s)ds)Q{x)^-[j {af{s) + b'{s))ds)&{x), 
we have proved the reconstruction formula (|4.ip . This finishes the proof of Thcorcin l4.2l (} 

Proof of Theorem \4-l\ We prove here that the parameters M, Q and A arc uniquely determined from 
the knowledge of the high energies of the transmission operator Tr. Note that the proof with the 
high energies of Tl is the same. Consider Tr^i and Tf{_2 two transmission operators corresponding to 
parameters Mj,Qj,Kj with j ~ 1,2 where moreover m,q ^ are supposed to be known and fixed. In 
what follows, we shall denote all the objects associated to T/j by the usual notations with a lower index 
j- 

We assume that Tr^i = Tr,2- From the definition of Fi{X) it follows then that -Fi,i(A) = Fi^2W- We 
identify now the terms of same orders in the asymptotic expansion (|4.ip . Since ipj 4* s-re dense in ?i, we 
get 

ei(a;) = 62(2;), VxGM (4.22) 
Ai{x) = A2{x), VxeM. (4.23) 

Let us analyse the term of order first. From (|4.22p and (|4.3p . we have 

— iPi — «(c+4 — 00,1)2: = —i(i2 — i{c+a ~ 00,2)2; + 2fc7r, Vx G M, (4.24) 

where k E Z. If we derivate (|4.24|) with respect to x, we thus obtain 

co,i - c+.i = Co, 2 - c+,2- (4-25) 

Hence using (|4.25|) and (|2.29p . we see that the quantity 

X = co-c+ = qQ^^±^, (4.26) 

is uniquely determined. 

We analyse now the term of order 0(A~^). From (|4.23p . (|4.3p and (|4.22p again, we have 

+00 /'+00 

W^{s)ds = / Wi{s)ds. (4.27) 

-00 J —00 

Using that (x) — of {x) + (x) and the expressions of the potentials ai and b given by (|2.15p and the 
definition of the Regge- Wheeler variable (|2.6p . we can compute explicitely the integrals that appear in 
(|4?27)) . In fact we have 

/ W^{s)ds = {l + -r( )+m'{r+-ro). (4.28) 

J-00 2 Vro r+/ 

By homogeneity in I and since m is considered as known and fixed, we deduce from (|4.27p and (|4.28p 
that 

= (4.29) 

ro,ir+^i fo,2?'+,2 

?'+.i-''o.i = r+,2-ro,2. (4.30) 
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Hence the quantities 



Y=''-^^, Z = r+-r^, (4.31) 



are uniquely determined. 

We can now show the uniqueness of the parameters M, Q and A as follows. We first note the following 
relation 

X = qQY. (4.32) 

Since X, Y are uniquely determined and q is supposed to be known and fixed, we deduce from (|4.32p that 
Q is uniquely determined, i.e. Qi = Q2 = Q- 

Moreover, from (|4.3ip we deduce that — tq and r^r^ are uniquely determined. Hence so are tq and 
r+ as the unique solutions of the obvious polynomial of second order. Now recall tq and are roots of 
F{r) = 0. The equations -F(?'o) = and F{r^) — can be written using (|2.2p as the linear system 



M 



The determinant of (|4.33p is | "^ror^ ^'^'^ clearly nonzero. Hence (M, A) are the unique solutions of the 
system (|4.33p whose coefficients depend only on tq , ?"+ , Q which are uniquely determined by the previous 
discussion. We thus conclude that M and A are also uniquely determined, i.e. Mi = AI2 and Ai = A2 
and the proof of Theorem 14. II is finished. <^ 



4.2 The inverse problem on an interval of energy 

In this last subsection, we solve the inverse problem when the reflection operators L or R are supposed 
to be known on a (possibly small) interval of energy. We follow the usual stationary approach of inverse 
scattering on the line and refer to [SI [6] for a presentation of the general method in the case of one- 
dimensonal Schrodinger operators and to |T] for an application to massless Dirac operators (see also 
|12[ I15j for massive Dirac operators) . We first determine a stationary representation of the scattering 
operator S expressed in terms of the usual transmission and reflection coefficients (here matrices) . We do 
this by a serie of simplications of our model which finally reduces to the exact framework studied in [1] . 
We then use the exponential decay of the potentials to show that the reflection coefficients R and L can 
be extended analytically to a small strip around the real axis. In consequence, the reflection coefficients R 
or L are uniquely determined on K if they are known on any interval of energy by analytic continuation. 
At last, we use the results of [1], a classical Marchenko method, to prove that the parameters M, Q and 
A are uniquely determined by the knowledge of i?(^) or i(^) for all energies. 
Recall that the scattering operator S is defined by 

S = {w+)*w-, 

where the global wave operators are given when A > by 

w^ = K^) + K+^r (4-34) 

with 

Wf" ,=s- lim e''"e-'*""P^, lim e'*" 6''*"+ P±. (4.35) 
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We now use the unitary transform U introduced in (|3.10ip and the corresponding simphfied expressions 
of obtained in (pHO]) and (jiJBll to express (|4.34p as 

= UW^{A, Ao) (^e'^''""=P^ + e*r'^e*^''=+^P±) . (4.36) 



Here we have used the notations introduced in Subsections l3.2l and l4.1l Let us denote by G± the operators 
gir caxp^ +6*^" ^e^^ appearing in (14.361) and by S{A, Aq) the scattering operator associated to the 

operators A and A^, i.e. 

S{A,Ao) = {W+{A,Ao)rw-{A,Ao). 
Using the unitarity of U we thus immediately get the foUowing expression for the scattering operator S 

S = GXS{A,Aq)G-. (4.37) 

The couple of operators {A,Aq) acting on Ti, turns out to fit the framework studied in [1]. Recall that 
they are given by = ^^D^ and A = Aq + W{x) where the potential W{x) = e*^''^-(^)(a(a;)r2 + 
6(x)r°)e~*'" C-(x) ^YiQ 4x matrix-valued function 



W{x) 



k{x) 
k*{x) 



Here k*{x) denotes the transpose conjugate of the matrix- valued function k{x). Moreover W satisfies 
()4.12|) and ()4.13|) and thus its entries belong to L^(M). This is precisely the kind of operators studied 
in Ij. Note however that our potential W is better than L-'^(]R) since it is exponentially decreasing at 
both ends x ztoo. This will be used hereafter. As a consequence, we can use the following stationary 
representation of S{A, Aq) obtained in [T]. Let us introduce the unitary transform J- on Ti defined by 

.FV(0 = ^ / e-'i^'-«^(a:)dx, (4.39) 
V 27r jR 

then we have (see [T], p 143)) 

SiA,Ao)=T*SoiO^, (4.40) 
where the scattering matrix 5*0 (^) takes the form 

Here Ti(^) and Tji{(^) are 2x2 matrix- valued functions which correspond to the usual transmission 
coefficients of S whereas L(^) and R{(,) are 2x2 matrix-valued functions which correspond to the usual 
reflection coefficients of S. We refer to Sections 2 and 3 of [1 for the definition and the construction of 
the scattering matrix 5*0 (0- Hence (|4.37p becomes 

S^iTG+rSoiO^G-. (4.42) 

We now finish our factorization of the scattering operator S as follows. Using 2x2 block matrix notations, 
we note that 

_ , e^f^ \ f 6*^=+^ \ _ f 1 \ f e^^o^i; q 
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and we define two unitary transforms F± on Ti. by 



and 



=r[ n .-"c+. ) m ^-^J{ g„c-.c+. ) ^{x)dx. (4.44) 



Then we have 



= ( M F+, = ( J 1 F_. (4.45) 



Hence we conclude from (|4.45p that the scattering operator (|4.42p factorizes as 



We summarize this result as a Proposition 

Proposition 4.1 The scattering operator S has the following stationary representation. If F± are the 
unitary transforms defined in ^-43^ o-nd ^.44^ , then 

S^F^S{^)F_, (4.47) 

where the 4x4 scattering matrix S'(^) is given by 

S{0 - y ^(^) e-^PTniO ) ' ^^■^^> 

and the quantities Tl,Tr and L, R are the 2x2 matrices that correspond to the transmission and reflection 
matrices of S{A, Aq) respectively and are obtained in ^Tj, Sections 2 and 3. 

Remark 4.1 As the notations suggest, the diagonal elements of the scattering matrix S{£^) given in 
{4-4^^ O'T^ simply the stationary representations of the transmission operators and T/j introduced in 
Section\^ 112. 33\) . The anti-diagonal elements of S{^) are in turn the stationary representations of the 
reflection operators L and R in \2.34^ . 

Remark 4.2 The unitary operators F± appearing in the stationary representation ^.47\ l of S are natural 
in the following sense. Let us define the two selfadjoint operators on Ti. 

H+ - {T^D^ + c+)P+ + {T^D, + co)P-, H- - {T^ + c^)P+ + {T^ + c+)P^. 

Hence it is clear from \4-34^ and {4-3^ that the global wave operators can be written in a classical form 
as 

W^ = s- lim e'*^e-'*^^ 

t — *±oo 

Now it is an easy calculation to show that the unitary transforms F± introduced in Il4-43\ l and \4-44^ f^f^ 
precisely the unitary transforms which diagonalize the operators , i.e 

= FIM^F±, 

where denotes the multiplication operator by ^. We conclude that J^.-^Tp together with ^4-4^^ 
expected stationary representation of the scattering operator S. 
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In the sequel, we shall use the explicit link between our scattering matrix S{^) and the scattering 
matrix Sq {£,) thoroughly studied in [1] in order to solve the inverse problem. Let us first briefly summarize 
some of the main results obtained in [T]. Under the assumption W G the scattering matrix So{^) 

is continuous for ^ S M and tends to I4 when ^ ±00. It is also unitary for each ^ G R (see [1], Thm 
3.1 for a proof of these statements and for other properties on So{^))). Moreover, the following partial 
characterization result holds: 

Theorem 4.3 (pLj, Thm 6.3) Assume that the reflection operators R{(,) and L{S) he 2x2 matrix valued 
functions satisfying 

sup ||i?(0|| < 1, sup ||L(0|| < 1, |l^(a)|l e Li(R), ||L(a)|| e L\R), (4.49) 

/'+00 />0 

/ a\\R{a)\\'^da< 00, / a\\L{a)\\'^ da < 00, (4.50) 

Jo J -00 

where R{a) and L{a) denote the usual Fourier transform of R{^) and L{^) and \\.\\ is the euclidean norm 
of a given matrix. Then the matrix-valued function k{x) G iy"'^(M) in ^.38^ ) (and thus the potential W{x) ) 
can he uniquely recovered from the knowledge of R{^) and L{^) for all ^ G M. 

We make several comments on this result and how we can apply it to our model: 

• The proof of the above theorem uses a classical Marchenko method. For instance, the matrix-valued 
function k{x) can be obtained after solving the following Marchenko integral equations for a > 
(see [J, eq. (6.9) and (6.11)) 

Bi{x,a) ^-R{a + 2x)+ / Bi{x,j)R{S + + 2xy R{a + 6 + 2x)d-yd6, (4.51) 

Jo Jo 

B2{x,a) = -L{a-2x)* + I / B2{x,-i)L{5 + - 2x)L{a + 5 - 2x)* d-^dS. (4.52) 

Jo Jo 

Under the assumption (|4.49p . the integral equations (|4.5ip and (|4.52p are uniquely solvable in 
L^(R+) ([1], Thm 6.2). Moreover, under the additionnal assumption (|4.50p . the matrix-valued 
function k{x) defined using the boundary values of i3i and B2 by the formulae (see [T], eq. (4.19)) 

k{x) = 2iBi{x, 0+), Vx > 0, k{x) = -2iB2{x, 0+), Vx < 0, 

can be shown to be in i^(R) and thus corresponds to the potential we are looking for. 



• If the potential W belongs to L^(R), then the condition (|4.49p is automatically satisfied (see [T], 
Thm 4.2 and eq. (6.17)). Although this condition is the natural one under which one could expect 
to reconstruct the potential k in the class L^, the authors of [1] had to add the extra assumption 
()4.50p (which must then be checked) in order to prove their result. We refer to [1], p. 154 for 
more details on this point. In our case, we shall prove the condition (|4.50|1 as follows. Using the 
exponential decay of W, we are first able to show that the reflection coefficients i?(^) and L(f ) (in 
fact the whole scattering matrix Sq{^)) are analytic on a small strip around the real axis. Moreover 
the functions R{. + irj) and L{. + irj) can be shown to belong to L^(R) uniformly for each \r]\ small 
enough. It follows then from standard results on the Fourier transform (see for instance p4 , Thm 
IX. 13) that R{a) and L{a) satisfy 

e'l"li?(a) G i2(R), e'l"lL(a) G L^(R), Ve smaU enough. 
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from which (|4.50p foUows immediately. 



• From (|4.5ip and (|4.52[) and the reconstruction procedure explained above, we see that the knowledge 
of R{£,) and L(^) for all ^ e M is used to recover the potential k{x) for all a; e R. In fact it is only 
enough to know either R{^) or L(^) for all ^ e M since then the whole scattering matrix 5o(C) can 
be uniquely recovered. The procedure is explained in [T], p. 147, eq.(5.3)-(5.5) and we reproduce it 
for completeness. Assume for instance that R{(,) is known for all ^ G M. Then the transmission 
coefficients T^i^) and Tj^,{S,) can be obtained performing the factorizations 

momo* ^h-R{omr, Tn{o*TR{o = h-R{o*m): ee^- (4.53) 

Under the assumption k e i^(M), it was shown in [l] that the above factorization problems are in 
fact left or right canonical Wiener-Hopf factorization in the Wiener algebra and thus lead to 
unique Tl(^) and Tji{£_) (see for instance [H], Thm 9.2, p. 831). At last, the reflection coefficient 
L{^) is recovered from by the formula 

m = -TRiom*iTL{o*r'- (4.54) 



• Eventually we explain how we can apply this result to our model. From Proposition [4711 we assume 
for instance that e~^'^^R{£_) is known for all ^ e M. Then it is easy to see from (|4.53|) and (|4.54p 
that we can uniquely recover Tl(^) and Tr{^) by performing Wiener-Hopf factorizations and then 
e^^^L{£) for all ^ € M. Note that the exponential term e^^*'^ disappears in the factorization (|4.53[) . 
If we assume that the assumptions ()4.49p and (I4.50p hold (this will be checked below) , then we can 
apply Thm l4.3l as follows. Multiplying the integral equations (|4.5ip and ()4.52|) by e~^*'^ and solving 
them, we conclude that we can uniquely recover e^*'^fc(x) (and not k{x)) for all a; G K. We shall 
show below that this implies the uniqueness of the parameters Af, Q and A of the black hole. 

Let us now show the analyticity of i?(^) and L{^) on a small strip around the real axis and prove 
there the uniform estimates mentioned above. To do this we need to introduce some objects whose 
existence has been shown in 1 , Sections 1, 2 and 3. The reflection coefficients R{C) and L(^) can be 
expressed in terms of solutions of the stationary problem 

T^D^ + W{x)\x{x,S,)=iX{x,£,), CeM (4.55) 

where A'(x,^) is understood as 4 x 4 matrix- valued function. Of special interest are the Jost solutions 
Fi(x,£,) and Fr{x,£,) of (j4.55p which are singled out by the specific asymptotics at infinity 

Fi{x,0^e'^'i^h + o{l)), x-^+oo, 

F,(x,C) =e'^'«"(/4 + o(l)), x^-oo. 

For each ^ £ K, these two solutions exist, are fundamental matrices of (|4.55p and are related as follows 
([1], Proposition 2.2). There exist two 4x4 matrix valued functions ai{£^) and ar{^) such that 

Fi{x,O^Fr{x,OMO, Fr{x,O^Fiix,0<^r{0, 

and satisfying ai{^)ar{£,) = CLr{£,)ai{£,) — I4 for all ^ e R. Note that Fi{x,S^) and Fr{x,^) satisfy the 
asymptotics (in the opposite ends) 

Fi{x,0=e'^'^''{ai{0+o{l)), x ^ -00, (4.56) 
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Fr{x,0 = e'^'^^ariO + X ^ +00. 

Let us now express ai(^) and ar{^) using 2x2 block matrix notations as 



Then the reflection coefficients are defined by ([T', eq. (3.6) and (3.7)) 

Since the situations are obviously symmetric, we shall only prove the analyticity and the uniform 
estimate on a small strip around the real axis for -R(^) (the proof for L(^) being identical). Moreover, 
we shall only consider the definition -R(^) = — aii(C)~^^^i2(C) for simplicity. To go further, we use some 
integral representations of the coefficients aii(^) and ai2{C) obtained in yy. These are given in terms of 
the Faddeev matrix Mi{x,0 defined by 

It is easy to see from (|4.55p that M/(x,^) must satisfy the integral equation (P, eq. (2.12)) 

p-\-oo 

Mi(x,i) iV^ / e-''i^'«(^--)W^(y)Mi(2/,C)e^^'«(^-"^dy, (4.57) 

J X 



and from (|4.56p that Mi{x,£^) must satisfy the asymptotics AIi{x,£,) =14 + o(l) when x — > +00. In fact, 
using once again 2x2 block matrix notations for AIi{x,S,) 



Mi{x,0 



Mn{x,0 Mi2{x,0 
Mi3ix,0 Mii{x,0 



and iterating (|4.57|) once, we get the uncoupled system of integral equations for AIi3{x,£^) and Mi4{x,£,) 
([T], eq. (2.15) and (2.16)) 



Mi3{x,0=i 



+ 00 



k{y)dy 



00 p-\-oo 



emy--)k{yrk{z)Mi3{z,0dzdy, 



k{y)*k{z)Mii(z, Cjdzdy, 



(4.58) 



(4.59) 



and similar equations for Mn(x,^) and Mi2{x,£) that we won't need. Eventually, the following integral 
representations for the coefficients a;i(^) and ai2{C) hold ([T], eq. (2.25) and (2.26)) 



aii{Cl = h-i / k{y)Mi3{y,0dy, 



0(2(0 = -i 



k{yyMu{y,Ody. 



(4.60) 



(4.61) 



We first study the coefficient a;2(C) expressed in terms of the Faddeev matrix Mi4{x,£,). Under the 
assumption k g L^(R), a solution Mi4{x,£,) of (|4.59p with the right asymptotics is easily shown to exist 
by iteration. Moreover for each fixed a; e R, this solution can be extended to a continuous function in the 
variable ^ when Im^ < and analytic when Im^ < ([T, Proposition 2.3). We prove now the following 
result 
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Lemma 4.1 Define the function P{x,^) = JJ*""" e^'^™^"^'||fc(y)||(i?/. Then there exists k > Q small enough 
such that 

(i) For all ^ satisfying \ImS^\ < k and for all a; G M, the function P{x,S,) is uniformly bounded. 

(a) For each fixed X £ M, the Faddeev matrix AIi4^{x,S,) can be extended analytically to the strip \Im^\ < k. 

Moreover, for each such it satisfies the estimate 

\\Mu{x,0\\ < Ccosh(P(x,0). (4.62) 

(Hi) For each fixed x G M, the derivative Ml^{x,S,) of the Faddeev matrix w.r.t. the variable x can be 
extended analytically to the strip \Im^\ < k. Moreover, for each such ^, it satisfies the estimate 

\\MUix,m\ <CsmHP{x,0)- (4.63) 

Proof: The first assertion is a direct consequence of the definition of P{x, ^) and (|4.13p (take for instance 
K = ^ where a is the positive number that appears in (|4.13p ). Solving (|4.59p by iteration leads to set 
Mi4{x,£,) = I]r=o"»(2^'6 with uo{x,£,) = h and 

r-\-oo r-\-oo 

Un{x,0^ / e'^'i(-''y'>k{y)*k{z)ur,-i{z,Odzdy, Vn > 1. (4.64) 

Jx Jy 

By induction we get the estimates 

lkn(x,OII< ^(^;y\ VneN. (4.65) 

Together with (i), this entails the second assertion. To prove the third one, we consider the serie of 
derivatives ''^niXjO- From (|4.64|) . note that 

<(x,e) = - / e-2*«("-=")fc(a;)*fc(z)M„_i(z,^)dzd2/. 

By induction and using (|4.65p . we get the estimates ||m^(x,^)|| < C '^^^rfli)! for all n > 1 from which 
we deduce (iii). 





Corollary 4.1 Let k the positive number defined in Lemma \4-.1\ The coefficient a/2(C) analytic on the 
strip \ImS\ < K. Moreover, it satisfies there the estimate 

\\ai2m\=Om~'), Id -ex.. (4.66) 

Proof: The analyticity on the strip \Im£^\ < k follows directly from (|4.6ip and Lemma |4?T1 To prove the 
second assertion, we integrate by parts in (j4.6ip . For all ^ with j/m^ < k, we obtain 

aniO = J^e-^'^y {k'iy)Mi4y,0 + kiy)MUiy,0)dy. (4.67) 

Since k' also satisfies the estimate (|4.13p and using Lemma HTTl again, we conclude that ||a(2(C)ll ^ 



We now study the coefficient aa(0 expressed in terms of the Faddeev matrix Mi3{x, ^). Once again under 
the assumption k S L"'^(R), a solution Mi3{x,£,) of (|4.58p with the right asymptotics is easily shown to 
exist by iteration. Moreover for each fixed a; G K, this solution can be extended to a continuous function 
in the variable ^ when Imf > and analytic when Im^ > (p]. Proposition 2.3). Using the same function 
P{x,^) and positive number k as in Lemma |4.1[ let us prove the following result 
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Lemma 4.2 For each fixed x € R, the Faddeev matrix Mi3{x,^) can be extended analytically to the strip 
|/to^| < K. Moreover, for each such ^, it satisfies the estimates 

\\Mi3{x,m < Ce2|^™«ll-lsinh(P(a;,0). (4.68) 
l!Mn(a:,OII <|J(l + e'"™«""'), iei>l. (4.69) 
Proof: We solve (|4.58p by iteration. Hence we set Mi3{x,£^) — J2'^=o "^"ni^ ' with 

r+oo 

vo{x,0^^ e2^«(^-^)fc(y)dy, 

and 



oo r + oc 



vn{x,0^ / e^'^^y-^h{yrk{z)vnMz,Odzdy. (4.70) 

We can prove the following estimate by induction 

which implies immediately (|4.68p . Moreover, since P{x, ^) is uniformly bounded on |/m^| < k, we deduce 
from (|4.68[) the analyticity of Mi3(x, £,) on the same strip. To prove (|4.69p . we integrate by parts in (|4.58p 
w.r.t the variable y. For all ^ with |/to^| < k, we obtain 

Mi,{x,0 = --^--^j^ e'^iy{k*Y{y)dy 



k*{x)K{x) e 



-2i£,x p+oo 



2ii 



e''^y[{k*y{y)K{y) ~ k* {y)k{y)Mi:,{y,o)dy, (4.72) 



where we have introduced the function K{x) = J^°° k{y)Mi3{y,£^)dy. Now using (|4.13p for k and fc', 
(j4.68p and the uniform estimate ||_ft'(a;)|| < C for all ^ with \Im^\ < k, we deduce that (|4.69p holds when 
|<^| is large from 





Corollary 4.2 Let k be the positive number defined in Lemma \4-.1\ Then the coefficient a/i(^) is analytic 
on the strip < k and tends to I2 when |^| —^ 00. Furthermore, possibly considering smaller k , the 

coefficient a/i(^) is invertible on the strip < k and aYi^{£_) is analytic and uniformly bounded there. 

Proof: The first assertion is a direct consequence of ()4.60p and Lemma l4?2] Since a;2(C) tends to I2 when 
1^1 00, a;2(0 is clearly invertible for |^| large enough. Since ai2{£,) is also invertible on the real axis 
([1], Proposition 2.10), we conclude that a;2(0 is invertible on a strip \Im^\ < e with < e < k small 
enough and that af^^{(,) is analytic and uniformly bounded on |/m^| < e. Denoting this e by k, we have 
proved the corollary. 



Let us put all these results together. Since i?(f) = —CLjl^{£^)ai2{^), Corollaries 14.11 and 14.21 implv that 
the reflection coefficient R{^) is analytic on a strip |/m^| < k where k is a small enough positive number. 
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Moreover, using the estimates of the same corollaries, we see that \\R{. + i7])\\ e L'^iM.) for aU |?7| < k. In 
fact, we have 

sup \\R{. + iri)\\L2 < oo. 
\v\<K 

Finally it follows from Thm IX. 13 in 24J that the Fourier transform R(a) satisfies the estimate 

e''l"l||i?(a)|| e L2(M). (4.73) 

In particular, the assumption ()4.50|) in Thm 14.31 is satisfied by 
We finish this paper solving the inverse problem. 

Theorem 4.4 Assume that one of the reflection matrices L{S^) or e~^*^i?(^) appearing in J^.^^p is known 
on a (possibly small) interval o/M. Assume moreover that the mass m and the charge q^Q of the Dirac 
fields are known and fixed. Then the parameters M, Q and A of the dS-RN black hole are uniquely 
determined. 

Proof: We only give the proof when the reflection matrix e~^*^i?(^) is supposed to be known on an 
interval / of M since the proof with L(^) can be treated the same way. We consider thus e~^''^^i?i(^) 
and e~^''^^i?2(0 two reflection matrices corresponding to parameters Mj, Qj and Aj with j — 1,2 where 
moreover the parameters m, g ^ are supposed to be known and fixed. As usual we shall denote all the 
objects related to e'^"^^' Rj{S,) by a lower index j in what follows. Assume that e~^*''ii?i(^) = e~^'^^^ R2{£,) 
for all ^ e /. By analyticity, we thus have 

Using the procedure explained after Thm 14.3] this also entails that 

Thanks to (|4.73p and the corresponding result for L(^), we can apply Thm l4.3l fand the remarks following 
this Theorem). Hence we obtain the equality e^'^^fci(a;) = e^'^^^k2{x) for all x e R or equivalently 

e^^^'P^Wi{x)^e^'^'^^W2{x), Vx e R. (4.74) 
Now recall that is a positive function since 

WHx)=ai{x)+b'{x) = (z + i)':^+m2F(r), 

Hence taking the square of ()4.74|) and then the modulus, we have 

Wl{x)^al^{x)+b\{x)^al^{x)+bl{x)^W^{x), Vx G R. (4.75) 
Note in particular that 

/ + 00 /• + 00 

Wl{s)ds = / Wi{s)ds. (4.76) 
-oo J ~oo 

Moreover by homogeneity in / and since ai and b are positive functions, we deduce from (|4.75p that 

ai,i{x)^ai^2{x), bi{x)^b2{x), Vx £ M. (4.77) 
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Now since 

by (|TTB)) it follows from and (|iT7)) that 

or equivalently that 

(ii - Cf (a;) = /32 - (a;) + fcyr, Vx e M, (4.78) 
where k E Z. Derivating ()4.78p . we obtain 

ci{x) = C2{x), Va;eR, (4.79) 

If we let tend x to ±00, we obtain from (|4.79p and (|2.15p 

co,i = Co, 2, c+,1 = (4.80) 

We notice eventually that (|4.76p and (|4.80p are precisely the conditions under which the parameters 
M, Q and A were shown to be uniquely determined in the proof of Theorem 14. II (see precisely the condi- 
tions (|4.25p and (|4.27p ). We thus apply the same procedure as before to end up the proof of the Theorem. 
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